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Abstract 

We present exact calculations of the partition function of the g-state Potts model on (i) open, (ii) cyclic, 
and (iii) Mobius strips of the honeycomb (brick) lattice of width L y = 2 and arbitrarily great length. In the 
infinite-length limit the thermodynamic properties are discussed. The continuous locus of singularities of the 
free energy is determined in the q plane for fixed temperature and in the complex temperature plane for fixed 
q values. We also give exact calculations of the zero-temperature partition function (chromatic polynomial) 
and W(q), the exponent of the ground-state entropy, for the Potts antiferromagnet for honeycomb strips 
of type (iv) L y = 3, cyclic, (v) L y — 3, Mobius, (vi) L y = 4, cylindrical, and (vii) L y = 4, open. In the 
infinite-length limit we calculate W(q) and determine the continuous locus of points where it is nonanalytic. 
We show that our exact calculation of the entropy for the L y = 4 strip with cylindrical boundary conditions 
provides an extremely accurate approximation, to a few parts in 10 5 for moderate q values, to the entropy 
for the full 2D honeycomb lattice (where the latter is determined by Monte Carlo measurements since no 
exact analytic form is known). 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase transitions and critical 
phenomena fjl |2j. On a lattice, or, more generally, on a (connected) graph G, at temperature T, this model 
is defined by the partition function 

{*»} 

with the (zero-field) Hamiltonian 

(y) 

where <Xj = 1, ...,q are the spin variables on each vertex (site) i € G; (3 = (ksT)^ 1 ; and (ij) denotes pairs 
of adjacent vertices. The graph G = G(V,E) is defined by its vertex set V and its edge set E; we denote 
the number of vertices of G as n — n{G) = \V\ and the number of edges of G as e(G) = \E\. We use the 
notation 

K = /3J , a = u- 1 =e K , v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., u > corresponding to oo > T > for the Potts ferromagnet, 
and (ii) < a < 1, i.e., — 1 < t; < 0, corresponding to < T < oo for the Potts antifcrromagnct. One defines 
the (reduced) free energy per site / = —@F, where F is the actual free energy, via 

f({G},q,v) = lim k[Z(G,«,«) 1 /«] . (1.4) 

n — ^oo 

where we use the symbol {G} to denote linin^oo G for a given family of graphs. In the present context, this 
n — > oo limit corresponds to the limit of infinite length for a strip graph of the honeycomb lattice of fixed 
width and some prescribed boundary conditions. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex set V and an edge 
set E' C E. Then Z(G, q, v) can be written as the sum 

z(G, ff ,«) = £ i k[G ' ) v <G ' ) 

G'CG 
n(G) e(G) 

= E E^g r w s (1.5) 

r=fe(G) s=0 

where fc(G') denotes the number of connected components of G' and z rs > 0. Since we only consider 



connected graphs G, we have k{G) = 1. The formula (1.5) enables one to generalize q from Z + to R+ 
(keeping u in its physical range). This generalization is sometimes denoted the random cluster model |3| ||]; 
here we shall use the term "Potts model" to include both positive integral q as in the original formulation 



in eqs. (1.1) and (1.2), and the generalization to real (or complex) q, via eq. (\L.aj) . The formula (1.5) shows 
that Z(G, q, v) is a polynomial in q and v (equivalently, a) with maximum and minimum degrees indicated in 
cq. (p..5j). The Potts model partition function on a graph G is essentially equivalent to the Tutte polynomial 
0-0 and Whitney rank polynomial 10 -j^] for this graph, as discussed in the appendix. 

In this paper we shall present exact calculations of the Potts model partition function for strips of the 
honeycomb (equivalently, brick) lattice of width L y = 2 vertices and arbitrary length equal to m bricks, with 



1 



boundary conditions of the following types: (i) (FBC y , FBC X ) — free or open, (ii) (FBC y , PBC X ) = cyclic, 
and (iii) (FBC y: TPBC X ) — Mobius, where FBC and (T)PBC refer to free and (twisted) periodic boundary 
conditions, respectively, and the longitudinal and transverse directions are taken to be x (horizontal) and y 
(vertical). This is a natural extension to the honeycomb lattice of our previous exact calculations of Potts 
model partition functions on finite- width, arbitrary-length strips of other lattices |l^]-]l7|. Note that while 
the Potts model partition functions, or equivalently, Tutte polynomials for the limit of infinite 2D triangular 
and honeycomb lattices are simply related by duality, this is not true for finite-width strips of the triangular 
and honeycomb lattices, owing to boundary effects. (Dual graphs are discussed further in the appendix.) In 
addition to these results for arbitrary temperature, we present exact calculations for the interesting case of 
the zero-temperature antiferromagnet for strips of the honeycomb lattice of the following types: (iv) L y = 3, 
cyclic, (v) L y = 3, Mobius, (vi) L y — 4, {PBC yi FBC x ) ^cylindrical, and (vii) L y — 4, free. 

There are several motivations for these exact calculations of Potts model partition functions for strips 
of the honeycomb lattice. Clearly, new exact calculations of Potts model partition functions are of value in 
their own right. From these, in the limit of infinite-length strips, one derives exact thermodynamic functions 
and can make rigorous comparisons of various properties among strips of different lattices, e.g., square, 
honeycomb, and triangular. In particular, with these exact results one can study both the T = critical 
point of the g-state Potts ferromagnet and T = properties of the Potts antiferromagnet. Further, we shall 
use our exact results to study how, for a given value of g, the ground state degeneracy per site (exponent of the 
ground state entropy per site) for the Potts antiferromagnet approaches its value for the infinite honeycomb 
lattice. We shall show that the exactly determined function for the ground state degeneracy calculated 
for the strip with an infinite-length strip with a rather modest width of L y = 4 and cylindrical boundary 
conditions is already an extremely good approximation, for a wide range of q values, to the corresponding 
function for the full 2D honeycomb lattice (the latter being determined by Monte Carlo measurements). This 
is very useful since no exact expression for the ground state degeneracy per site for q > 2 is known for the 2D 
honeycomb lattice. In addition, it was shown |]l4| that these calculations can give insight into the complex- 
temperature phase diagram of the 2D Potts model on the given lattice, which is again useful since, except for 
the q = 2 Ising case, the latter phase diagrams are not known exactly. Finally, via the correspondence with 
the Tutte polynomial, our calculations yield several quantities of relevance to mathematical graph theory. 

Various special cases of the Potts model partition function are of interest. One special case is the zero- 
temperature limit of the Potts antiferromagnet. For sufficiently large g, on a given lattice or graph G, this 
exhibits nonzero ground state entropy (without frustration). This is of interest as an exception to the third 
law of thermodynamics jl^, ^9|. This is equivalent to a ground state degeneracy per site (vertex), W > 1, 
since So — kslnW. The T — (i.e., v = —1) partition function of the above-mentioned g-state Potts 
antiferromagnet on G satisfies 



where P{G, q) is the chromatic polynomial (in g) expressing the number of ways of coloring the vertices 
of the graph G with q colors such that no two adjacent vertices have the same color || [l(], ||(| pifl . The 
minimum number of colors necessary for this coloring is the chromatic number of G, denoted x(@)- We have 



Z(G,q,-l) = P(G,q) 



(1.6) 



W({G},q)= lim P(Cg) 1 / 



(1.7) 



At certain special points q s (typically g s = 0, 1, ..,x(G)), one has the noncommutativity of limits 



lim lim P(G,g) 1/n ^ lim limP(G,g) 1/ 



(1.8) 
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and hence it is necessary to specify the order of the limits in the definition of V7({G}, q s ) S 



Using the formula (L5) for Z(G, q, v), one can generalize q from Z + not just to R+ but to C and a from 
its physical ferromagnetic and antiferromagnetic ranges 1 < a < oo and < a < 1 to a £ C. A subset of 
the zeros of Z in the two-complex dimensional space C 2 defined by the pair of variables (q, a) can form an 
accumulation set in the n — > oo limit, denoted B, which is the continuous locus of points where the free 
energy is nonanalytic. This locus is determined as the solution to a certain {G}-dependent equation jl3|. [l4|| . 
For a given value of a, one can consider this locus in the q plane, and we denote it as B q ({G}, a). In the 
special case a = (i.e., v = — 1) where the partition function is equal to the chromatic polynomial, the zeros 
in q are the chromatic zeros, and B q ({G},a = 0) is their continuous accumulation set in the n — ► oo limit. 
In previous papers we have given exact calculations of the chromatic polynomials and nonanalytic loci B q 



for various families of graphs; some of these are related earlier works are |22|]-||26|-|]48[ (see 14 for further 
refernences) . With the exact Potts partition function for arbitrary temperature, one can study B q for a ^ 
and, for a given value of q, one can study the continuous accumulation set of the zeros of Z(G, q, v) in the 
a plane; this will be denoted B a ({G}, q). It will often be convenient to consider the equivalent locus in the 
u = 1/a plane, namely B U ({G}, q). We shall sometimes write B q ({G}, a) simply as B q or B when {G} and a 
are clear from the context, and similarly with B a and B u . One gains a unified understanding of the separate 
loci B q {{G}) for fixed a and B a {{G}) for fixed q by relating these as different slices of the locus B in the C 2 
space defined by (q, a) as we shall do here. 

Following the notation in |^2f and our other earlier works on B q ({G}) for a = 0, we denote the max- 
imal region in the complex q plane to which one can analytically continue the function W({G},q) from 
physical values where there is nonzero ground state entropy as R\ . The maximal value of q where B q 
intersects the (positive) real axis was labelled q c ({G}). Thus, region Ri includes the positive real axis for 
q > q c ({G}). Correspondingly, in our works on complex-temperature properties of spin models, we had 
labelled the complex-temperature extension (CTE) of the physical paramagnetic phase as (CTE)PM, which 
will simply be denoted PM here, the extension being understood, and similarly with ferromagnetic (FM) and 
antiferromagnetic (AFM); other complex-temperature phases, having no overlap with any physical phase, 
were denoted Oj (for "other" ) , with j indexing the particular phase |49[ Q . Here we shall continue to use 
this notation for the respective slices of B in the q and a or a planes. Various special cases of Z(G,q,v) 



applicable for arbitrary graphs G were given in 1 14 



Just as one must take account of the noncommutativity in the definition of W, eq. (fL8|), so also one must 
take account of noncommutativity in the definition of the free energy: at certain special points q s (typically 
q s = 0, l...,x(G)) one has 

lim lim Z(G,q,v) 1/n £ lim lim Z(G, q, v) 1/n . (1.9) 

n — >oo q — >q s q — »g s n — >oo 



Because of this noncommutativity, the formal definition (1.4) is, in general, insufficient to define the free 
energy / at these special points q s ; it is necessary to specify the order of the limits that one uses in eq. ([O]) . 
We denote the two definitions using different orders of limits as f qn and f nq : 

f m {{G},q,v)= lim lim n' 1 \n Z(G, q,v) (1.10) 

n — >oc q^-q s 

fqn({G},q,v) = lim lim n' 1 In Z(G, q, v) . (1.11) 

q — >q s n — >oo 



In Ref. ]22[ and our subsequent works on chromatic polynomials and the above-mentioned zero-temperature 



antiferromagnetic limit, it was convenient to use the ordering W({G},q s ) — lim 9 ^ gs lim^oo P{G, q) 1 ^ 
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since this avoids certain discontinuities in W that would be present with the opposite order of limits. In 
the present work on the full temperature-dependent Potts model partition function, we shall consider both 
orders of limits and comment on the differences where appropriate. Of course in discussions of the usual 
g-state Potts model (with positive integer g), one automatically uses the definition in eq. (1.1) with (1.2) 
and no issue of orders of limits arises, as it does in the Potts model with real q. Consequences of the 
noncommutativity (1.9) have been discussed before jlj, [l(| [l7| . 

As derived in [|l4| , a general form for the Potts model partition function for the strip graphs G s considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated subunits (e.g. the 
columns of squares of height L y vertices that are repeated L x — m times to form an L x x L y strip of a 
regular lattice with some specified boundary conditions), is 

Nz,G s ,X 

Z((G s ) m ,q,v)= ry,; ,,(A Z , G j(q,v)) m (1.12) 

3=1 

where Nz,g 3 ,\i cz.G s .ji an d ^z.G 3 .j depend on the type of recursive family G s (lattice structure and boundary 
conditions) but not on its length m. For special case of the T = antiferromagnet, the partition function, 
or equivalently, the chromatic polynomial P{{G s ) m ,q) has the corresponding form 

N P,G 3 ,X 

P((G s ) m ,q) = v,.,; A\r.<; .J'/.rii'" . (1.13) 

3=1 

For sufficiently large integer q values, the coefficients can be interpreted as multiplicities of eigenvalues of 
coloring matrices , and the sums of these coefficients are thus sums of dimensions of invariant subspaces 
of these matrices. For a given family G s of strip graphs we shall denote these sums as 

N Z,G 3 ,\ 
3=1 



and 



N 



P.G„ 



3=1 

Here we distinguish the quantities Nz,g s ,X> ^Vp,G s ,a; cz,g 3 j, cp,G a ,j', ^z,G a ,j, and Ap,G s ,j- Below, for brevity 
of notation, we shall sometimes suppress the Z or P subscripts and the s subscript on G s when the meaning 
is clear from context and no confusion will result. 

The Potts ferromagnet has a zero-temperature phase transition in the L x — > oo limit of the strip graphs 
considered here, and this has the consequence that for cyclic and Mobius boundary conditions, B passes 
through the T = point u = 0. It follows that B is noncompact in the a plane. Hence, it is usually more 
convenient to study the slice of B in the u = 1/a plane rather than the a plane. For the ferromagnet, since 
a — > co as T — > and Z diverges like 

a e(G s ) m 

this limit, we shall use the reduced partition function Z r 

defined by 

Z r ((G s ) rn , q, v) = a^ G ^Z((G s ) m , q, v) = u e « G ^ Z{{G s ) m , q, v) (1.16) 

which has the finite limit Z r — > q as T — > 0. For a general strip graph {G s ) m of type G s and length L x = m, 
we can write 

N Z,G 3 ,X N Z ,G 3 ,\ 

Z r ((G s ) m ,q,v) = M £ (( G =) m ) J2 </.<: .J'" £ rx.<; .,^z.<; .,.„;'" (1-17) 

3=1 3 = 1 
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with 

\z,G a , jt K=u e « G ^V m \ ZiGa!j . (1.18) 
For the Potts model partition functions of the L y = 2 cyclic and Mobius strips of the honeycomb lattice to 



be discussed below, the prefactor in eq. (1.18) is u 5 . 

For the following, it will be convenient to define some general functions. First, we define the following 
polynomial: 

^«>-^ -§-')■(*; 0^ 

and P{C n , q) is the chromatic polynomial for the circuit (cyclic) graph C n with n vertices, given by P(C n , q) 
(q-l) n + (q-l)(-l) n . 

Second, we define coefficients of degree d in q which are related to Chebyshev polynomials: 



(1.19) 



-"■ 2 J 

where U n (x) is the Chebyshev polynomial of the second kind, defined by 

71 ~ A O.n-X*-* 



U 2d {4) (1-20) 



* 3 



U n ( X )=J2(- 1 Y[ ./)(2x) n - 2j (1.21) 



where in eq. ( 1.21 ) the notation [S] in the upper limit on the summand means the integral part of ^. The 
first few of these coefficients are 

c (°) = l, c^=q-l, (1.22) 
c (2) = q 2 -3q+ 1, (1.23) 

and 

c ( 3 ) = q 3 ~5q 2 + 6q-l . (1.24) 



2 Chromatic Polynomial for the L y = 2 Cyclic and Mobius Strips 
of the Honeycomb Lattice 

Before giving our calculation of the Potts model partition functions for these strips, it is useful to discuss 
a particularly interesting special case, namely that of the zero-temperature antiferromagnet. For the open 
strip of m bricks, which will be denoted S m for short, 

P{hc, 2 X m, FBCy, FBC Xl q) = q{q - l){D 6 ) m = q{q - l){q 4 - 5q 3 + 10q 2 - lOq + 5) m (2.1) 

so that in the m — * oo limit, 

W{hc, 2 x oo, FBC yi FBC x ,q) = (q 4 ~ 5q 3 + lOq 2 - lOq + 5) 1/4 . (2.2) 

Since this has only discrete branch point singularities at the zeros of Dq, the continuous locus B = 0. 

The chromatic polynomials for width L y = 2 cyclic and Mobius strips of the honeycomb lattice can 
be obtained as special cases of the calculations in |}8| for homeomorphic expansions of cyclic and Mobius 
strips of the square lattice. These homeomorphic expansions involved the addition of k — 2 degree-2 vertices 



5 



on the upper and lower horizontal edges of the strips, thereby making cyclic or Mobius strips of p-gons 
with p = 2k such that two successive p-gons intersected on one common edge. Here, the degree A (i.e., 
coordination number) of a vertex in a graph G is defined as the number of edges (bonds) that connect to it. 
The case k = 3, i.e., p = 6, was thus the width 2 cyclic or Mobius strip of the honeycomb lattice, depending 
on the longitudinal boundary conditions. We denote these strips as L m = he, 2 x m, FBC y , PBC X and 
ML m = he, 2 x m, FBC y , TPBC X for short. For quantities that are independent of the value of m in L m 
or ML m , we shall use the notation L and ML. The graphs L m and ML m have n — 4m vertices. For m 
large enough to avoid degenerate cases (i.e. to > 2 for the cyclic strip and to > 1 for the Mobius strip), the 
chromatic number for the L y x to strips of the honeycomb lattice with cyclic boundary conditions is 

x(hc,L y xm,FBCy,PBC x ) = 2 (2.3) 

while for the Mobius strips we have 

X (hc, L y x to, FBCy, TPBC X ) = 3 (2.4) 

independent of m in the indicated ranges. 

The results for the chromatic polynomials (in our current notation) are |3^| 

4 

P(hc, 2 x to, FBCy, PBC x ,q) = J2 cp,L J (Ap, Lj ) m (2.5) 

i=i 

and 

4 

P(hc, 2 x m, FBC y ,TPBC x ,q) = £ c P , M L,i(Ap > L, i ) m (2.6) 

3=1 

where, as indicated, the Ap,G,j's are the same for cyclic and Mobius strips |3£| 

Ap,l,i = 1 (2.7) 

Ap,l,2 = Di + D 3 = {q- l) 2 (2.8) 

\ P: L,3 = D 4 - D 3 = q 2 -4q + 5 (2.9) 

Ap,l,4 = D 6 = q 4 - 5q 3 + 10q 2 - lOq + 5 . (2.10) 
The corresponding coefficients are 

cp,l,\ - c& (2.11) 

CP.L.2 = CP,L, 3 = C (1) (2.12) 
CP,L,4 = - 1 (2.13) 

cp.mla = -c (0) = -1 (2.14) 

Cp.ML,2 = -C^ (2.15) 

c p ,ml,3 = c (1) (2.16) 

cp,mla = c (0) = 1 . (2.17) 

The equality of the coefficients (here, cp^L,i and cp^ML,i) °f the Xa,j that is dominant in region R\, inde- 
pendent of longitudinal boundary conditions, is a general result (i^, [T^| . 
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Fi gure 1: Locus Sq for the m — ► oo limit of the cyclic L y = 2 strip {-£/} of the honeycomb lattice with a = 0. Zeros of 
Z(L m ,q,v = —1) = P(L m ,q) for m = 20, i.e., n = 80 vertices, are shown for comparison. 

The singular locus B was given as Fig. 1 (a) in ]3c| . For our current generalization to nonzero temperature, 
it will be necessary to refer to details of this locus, so we show it as Fig. [y. As is evident in this figure, the 
locus (boundary) B separates the q plane into six regions: (i) R\ containing the real intervals q > q c = 2 and 
q < and extending outward to complex infinity, (ii) an innermost region i?2 containing the real interval 
< q < q c , and two complex conjugate pairs of regions, (hi) i?3, R% to the upper and lower right of q c , and 
(iv) i?4,i?! to the upper and lower left of q c . The dominant terms in these regions are (i) \p,l,a in Ri, (ii) 
Ap.L.3 in i?2, (hi) ^p,l,2 in i?3,i?3, and (iv) \px.i m Ri, Rl- The point q c is a multiple point where six 
curves, forming three branches of B, cross. This corresponds to the fact that at q c there is a degeneracy in 
magnitude of all of the Ap^j 's. The angles at which the various boundary curves emanates from q c are: 
(a) 7r/4 for the curve separating i?i and R$, (b) ir/2 for the curve separating i?3 from R4, and (c) 37r/4 for 
the curve separating R4 from i?2, and so forth for the complex-conjugate curves. The density of chromatic 
zeros is small on the boundary separating R2 and i?4, and its complex conjugate boundary. The locus B has 
support for Re(q) < as well as Re(q) > 0. Further, some chromatic zeros for hnite-m strips and a portion 
of the locus B for the m — > 00 limit both lie within the circle \q — 1| = 1. This is also true for the L y = 3 
cyclic and Mobius strips of the honeycomb lattice (see below). In this context, we recall that the conjecture 
later made in |4^] that chromatic zeros and, in the limit of infinite-length, the locus B, do not lie in the 
interior of the circle \q — 1| was stated to be specific to strips of the square and triangular lattice. 

The W function in these regions is given by 



W = (A0 1/4 for q € Ri 
\W\ = \q 2 -4 9 + 5| 1/4 for q £ R 2 



(2-18) 



(2.19) 
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|T^| = |<7-1| 1/2 for qeR 3 ,Rl (2.20) 

\W\ = 1 for q £ Ri,Rt . (2.21) 

Recall that for regions other than R\, only the magnitude of W can be determined unambiguously p2| . 
Note that W(q — 2) = 1. We next generalize this study to the case of finite temperature for both the 
antifcrromagnet and the ferromagnct. 

3 Potts Model Partition Function for Open L y — 2 Strip of the 
Honeycomb Lattice 

By using an iterative application of the deletion-contraction theorem for Tutte polynomials and converting 
the result to Z, or by using a transfer matrix method, one can calculate the partition function for the open, 
cyclic, and Mobius strips of the honeycomb lattice of width L y and arbitrary length m bricks, which we shall 
again denote as S m , L m , and ML m . We have used both methods. These calculations are a natural extension 
of our earlier calculations of Potts model partition functions for finite- width lattice strips of arbitrarily great 
length |14fl-||l7|. We begin with the open strip. The results are conveniently expressed in terms of a 
generating function 

oo 

T z (S m ,q,v) = Z(S m ,q 7 v)z m . (3.1) 



m=0 



We find 



where 



q, v) = — r (3.2 

(1 - A2,S,l2)(l - \ZS,2Z) 

A S fl=q(q + v) (3.3) 

A s ,i = -q 2 v 4 (l + v) (3.4) 

*z,s,(i,2) = l(T s ±^) (3.5) 

T s = v 5 + 6v 4 + lOqv 3 + 10q 2 v 2 + 5q 3 v + q 4 (3.6) 
R s = 4qv s - Aq 2 v 7 - 6q 3 v 6 - 2q 4 v 5 + v 10 + 32i; 8 + 8v 9 + 104w 7 g 

+ I96v 6 q 2 + 2Mv 5 q 3 + 208v 4 q 4 + q 8 + I0q 7 v + A5q 6 v 2 + 120q 5 v 3 . (3.7) 

The equivalent Tutte polynomial T(S m ,x,y) is given in the appendix. 

We have studied the locus B in the C 2 space defined by the variables (q 7 v) (or equivalently, (q 7 a)). 



with 



For a — 0, Z reduces to the chromatic polynomial (2.1), and the resultant W function, given in eq. (2.2) 



has only discrete branch point singularities at the zeros of Dq (namely, q ~ 0.690983 ± 0.9510565i and 
1.809017 ± 0.587785i), and consequently B = 0. As a initially increases from 0, B q consists of arcs that 
spread out from these four points. An illustrative case, a = 0.1, is shown in Fig. ^|. As a — * 1, these arcs 
shrink in toward the origin, q = 0. 

In Figs. H and || we show B u for q = 2 and 3, respectively. The positive real u axis, and the region 
reached by analytic continuation from it, form the (CTE)PM phase. Since the infinite-length, finite-width 



8 



\ 



0.5 



\ 



lm(q) 



-0.5 



-1 







0.5 



1.5 



2 



Re(q) 



Figure 2: Locus Bq for the n — * oo limit of the L y = 2 honeycomb strip {5} with free longitudinal boundary conditions, for 
a = v + 1 = 0.1. Zeros of Z(S m , q,v = —0.9) for m = 21, i.e., n = 86 vertices, are shown for comparison. 

strips are quasi-one-dimensional systems and the spin-spin couplings are short-range (specifically, nearest- 
neighbor), an elementary application of a Peierls argument shows that there is no long-range order and thus 
no broken-symmetry (FM or AFM) phases at any finite temperature. Related to this, the singular locus B u 
docs not cross the positive real u axis. 

4 Thermodynamic Properties of the Potts Model on the 2 x oo 
Strip of the Honeycomb Lattice 

We next discuss the physical thermodynamic properties of the g-state Potts model on the 2 x oo strip of 
the honeycomb lattice. These are independent of the longitudinal boundary conditions. The reduced free 
energy per site is given (with the f nq definition) by 



This is analytic for all finite temperature for both signs of the spin-spin coupling (ferromagnetic and anti- 
ferromagnetic). This analyticity property is equivalent to the above-mentioned fact that the singular locus 
B u does not cross the positive real u axis. It is straightforward to calculate from this the internal energy U 
and specific heat C per site. These have the high-temperature expansions 



/ = -lnA z ,s,i • 



(4.1) 




(4.2) 
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Figure 3: Locus B u for the m — > oo limit of the L y = 2 strip of the honeycomb lattice with free longitudinal boundary 
conditions, {S} with q = 2. Zeros of Z(S m ,q = 2,-u) in u for m = 21 (so that Z is a polynomial of degree 106 in u and hence, 
up to an overall factor, in u) are shown for comparison. 
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Figure 4: Locus B u '- same as in Fig. H for q = 3. 
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c = 



5k B (q-l)K 2 
iq 2 



1 + (£_?)# + 0^ 



(4.3) 




Figure 5: Specific heat C/ks for the g-state Potts ferromagnet on the 2 X oo strip of the honeycomb lattice as a function of 
K = J/(kgT). From bottom to top the curves are for q = 2,3,4. 

In Fig. H we show a plot of the specific heat for the ferromagnetic case. One sees that as q increases, the 
value of K at which the maximum in the specific heat occurs increases; i.e. this maximum occurs at lower 
temperature. This is a consequence of the fact that the spins are "floppier" for larger q, and hence one must 
cool the system to lower temperatures in order to achieve short-range ordering (as noted above, there is, of 
course, no long-range ordering). Further, as q increases, the height of the maximum increases. This can be 
interpreted as a result of the fact that as one makes q larger, one increases the effective number of degrees 
of freedom per site, and it is the onset of short-range ordering of these degrees of freedom that produces the 
maximum in the specific heat. 

It is also of interest to compare our exact calculations of the specific heat for the 2 x oo strip of the 
honeycomb lattice with our previous calculations for the 2 x oo strips of the square and triangular lattices. 
In Fig. H we show this comparison for the Ising case q = 2. One sees that the value of K at which the 
maximum occurs increases as one goes from triangular to square to honeycomb lattices. This is a consequence 
of the fact that the short-range ordering in the spins is strengthened as one increases the effective coordination 
number (vertex degree) of the lattice strip. The cyclic L y = 2 strips of the square and triangular lattices are 
A-regular, i.e., each vertex has the same degree, namely A = 3 and A = 4, respectively. This is also true of 
the infinite-length limit of the open L y = 2 strips of the square and triangular lattices. For the L y = 2 strip 
of the honeycomb lattice, and more generally, for a graph G and its n — > oo limit {G}, it is useful to define 
the effective coordination number (degree) as 

For the infinite-length limit of the L y = 2 strip of the honeycomb lattice, half of the vertices have degree 
2 and half have degree 3, so that the effective coordination number is 2.5. Thus, since one increases the 
(effective) coordination number from 2.5 to 3 to 4 going from the honeycomb to square to triangular L y = 2 
strips, one can achieve short-range order at a progressively higher temperature (smaller value of K). The 
same qualitative comparison applies for q = 3,4, etc. 
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Fi gure 6i Comparison of the specific heat C/kg (vertical axis) for the g-state Potts ferromagnet on the 2 X oo lattice strips 
for the Ising case, q = 2. In order of increasing position K m ax of the maxima in these curves, they apply for the (a) triangular, 
(b) square, and (c) honeycomb. 



Since the Potts ferromagnet has a T = critical point on infinite-length, finite-width strips, it is worth- 
while to examine the critical singularities. The specific heat has the low-temperature behavior 



C = 6k B K 2 (q-l)e 



-IK 



l + ^( q + 2)e- K + 0(e-™) 



-2K\ 



as K — > oo 



(4.5) 



Thus, as is typical of systems at their lower critical dimensionality and was true of the other finite-width, 
infinite- length strips for which we have performed exact calculations jlj, [if], |l7j , the specific heat exhibits an 
essential zero at T = as a function of temperature. In general, at a critical point the ratio p of the largest 
subdominant to the dominant Aj's determines the asymptotic decay of the connected spin-spin correlation 
function and hence the correlation length 



i = - 



i 

hip 



Since Xz.l,5 and \z,l,3 are the dominant and leading subdominant Az,Gj's, respectively, we have 

Az,L,3 



Pfm 



A 



Z.L,5 



(4.6) 



(4.7) 



and hence for the ferromagnetic zero-temperature critical point we find that the correlation length diverges, 
as T — > 0, as 

€fm - q^e 2K + 0(e K ) , as K -» oo . (4.8) 



Again, this is an essential singularity (divergence) as a function of temperature, as is generic for systems at 
their lower critical dimensionality. The detailed form of the exponent is similar to the behavior £ ~ q~ 1 e 2K 
of the correlation length at the T — critical point of the Potts ferromagnet on the L y = 2 infinite-length 
strip of the square lattice pfl . While the correlation lengths for other ID or quasi-lD systems also exhibit 
divergences with essential singularities, the detailed form of the exponent (coefficient of if) varies. Thus, for 
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example, £fm ~ q~ 1 e K as K — > oo for the ID Potts ferromagnet, while £fa/ ~ q~ 1 e 3K and £fm ~ q~ 1 e AK 
as if -> oo for the Potts ferromagnet on the L y = 2 strips of the triangular lattice JT(| and on the square 
lattice with next-nearest-neighbor spin-spin couplings, respectively. 

Except for the q = 2 Ising case, where, as is true on any bipartite lattice, there is an elementary 
equivalence between the ferromagnetic and antiferromagnetic signs of the spin-spin coupling J, the q-state 
Potts antiferromagnet on the 2 x oo strip of the honeycomb lattice is not critical at T — 0. This is reflected 
in the fact that B q only crosses the positive real axis at q = 2. For q ^ 2, the low-temperature behavior of 
the specific heat is given by C ~ const, x a as a — ► 0, while for q = 2, C ~ const, x a 2 as a — > 0, as follows 



from the equivalence and eq. (4.5). A plot of the specific heat for this antiferromagnetic case is given in Fig. 




2 4 -K 6 8 10 



Figure 7: Specific heat C/ks for the g-state Potts antiferromagnet on the 2 X oo strip of the honeycomb lattice as a function 
of — K = — J I {kgT) > 0. From top to bottom, the curves are for q = 2, 3, 4. 



5 Potts Model Partition Function for Cyclic and Mobius L y = 2 
Strips of the Honeycomb Lattice 

5.1 Results for Z 

For the cyclic and Mobius L y — 2 strips of the honeycomb lattice of length m bricks (G s ) m = L m , ML m , 
using the same methods as we did for the open strip discussed above, we calculate the partition functions 

6 

Z((G s ) m ,q,v) = J2 c z,GsA x z,GsA<l, v )) m fOT G S = L,ML (5.1.1) 

3=1 

where 

^Z,L,j = ^Z,ML,j , J = 1, 6 (5.1.2) 
A z ,l,i=« 4 (5.1.3) 
Az,L,2 =v 2 (q + v) 2 (5.1.4) 
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A 



Z,L,(3,4) 



V 

~2 



with 



<T + 4qw + Qv* + ± V#34 
i? 34 = 40<2« 3 + 24gV + 8q 3 v - 2q 2 v 3 + 8v 5 + 32v 4 + g 4 + v 6 



and 



For the cyclic strip the coefficients in eq. ( |5.1.1| ) are 



Cz,£,l = C (2) = q 2 - 3q + 1 



cz,Lj r = c (1) = q - 1 for j = 2, 3, 4 



.,(0) 



For the Mobius strip, 



czmla 



1 for j — 5, 6 



,(o) 



CZ,ML,2 = -C (1) 

Cz,ML,j = c (1) for J = 3,4 

CZ,ML,3 = 1 for J = 5, 6 . 



(5.1.5) 
(5.1.6) 
(5.1.7) 

(5.1.8) 

(5.1.9) 
(5.1.10) 

(5.1.11) 
(5.1.12) 
(5.1.13) 
(5.1.14) 



These are the same respective coefficient functions that we found for g-state Potts model partition function 
for the width L y = 2 cyclic and Mobius strips of the square lattice and for the cyclic strip of the triangular 
lattice (for the Mobius strip of the triangular lattice, the coefficients are algebraic rather than polynomial 
functions of q) O, n4 Hq] . The respective sums of the coefficients are 



and 



Cz,ml — q ■ 



(5.1.15) 



(5.1.16) 



In the T = (i.e., v = — 1) case for the Potts antiferromagnet, where the partition function reduces 
to the chromatic polynomial, two of the six \z,L,j's vanish and the others reduce to those for the L y = 2 
cyclic/Mobius strip given in pq] and above: 



<Vz,L,l = ^P,L,1 for v = -1 

Az,l,2 = Ap,l,2 = (q - I) 2 for v = -1 
\zx,3 = Ap,l,3 = q 2 ~ 4q + 5 for v = - 
^z.L.j =0 for v = — 1 and j = 4, 6 



Az,L,5 = A P , L ,4 = Aj = q ~ 5q 3 + 10q 2 - Wq + 5 for « = -1 
Here the sums of the corresponding coefficients are 

Cp,l = q(q - 1) 

and 



C 



P. ML 



= . 



(5.1.17) 

(5.1.18) 
(5.1.19) 
(5.1.20) 
(5.1.21) 

(5.1.22) 
(5.1.23) 
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For infinite temperature, v — 0, we have 

Az,L,5 = g 4 for v = (5.1.24) 

Xz,L,j=0 for v = and j = 1,2,3,4,6 (5.1.25) 
so that the partition function reduces to 

Z(L m , q, 0) = Z(ML m , q, 0) = (A^L, 5 ) m = q Am = q n ■ (5.1.26) 

From our exact calculations, it follows that in the full C 2 space spanned by (q,v), 

B({L})=B({ML}) . (5.1.27) 

This is the same result that we obtained for the analogous strips of other lattices |l^, [l4]. [Tg|. O] and is in 
accord with the conclusion that the singular locus is the same for an infinite-length finite-width strip graph 
for given transverse boundary conditions, independent of whether the longitudinal boundary condition is 
cyclic or Mobius. For B q we have proved this in by the use of crossing-subgraph strips, which subsume 
cyclic and Mobius strips as special cases (see further below). (Of course, the locus B is different if one uses 



open, as contrasted with cyclic/Mobius, longitudinal boundary conditions.) Owing to the equality (5.1.27), 
we shall henceforth, for brevity of notation, refer to both B({L}) and B({ML}) as B({L}) and similarly for 
specific points on B, such as q c ({L}) = q c ({M L}) , etc. 

5.2 B in the q Plane 

The locus B exhibits a number of qualitative features similar to those that we have discussed in earlier 
works [[l3| |lj, [u| |l7l] . Let us first consider the slice of this locus in the q plane, B q as a function of v, or 
equivalently, a = v + 1. As a increases from 0, the single crossing of B q which had been at q = 2 for a = 
splits into two crossings. At the larger one, the boundary is determined by the equation of the degeneracy of 
magnitudes of leading terms |Az,l,5| = IA.2.L.2I, which yields the value of q c as the real solution of the cubic 
equation q 3 + 3vq 2 +v 2 (2 — v)q + 2v 3 = 0. The second crossing is determined by the equation of degeneracy 
of magnitudes of | Ax,i,j | for j — 1, 2, 4, 6 and occurs, for < a < 1, at 

qinner = —2v = 2(1 - a) . (5.2.1) 

or equivalently, a = — {qinner — 2)/2. Parenthetically, we note that this is the value for q c , as a function of a, 
for the infinite- length limit of the circuit graph (5l], In the interval < q < qinner, ^z.l.3 is dominant, 
and, as was the case at a = 0, the crossing of B q at q = is determined by the equation of degeneracy of 
leading terms IA^l^I = IA^l^I. For sufficiently small a near 0, there is also a pair of complex-conjugate 
regions where \z,l,\ is dominant. The outermost portion of B q crosses the real axis at q — q c and q = 0; 
outside of this is the region including the real intervals q > q c and q < 0. As illustrations, we show B q 
in Figs. |s|-p~0| for a = 0.1, 0.5, and 0.9, respectively. For the case a = 0.9 we note that (a) the boundary 
includes cusp- like structures at q c and qinner, and (b) in the pair of small complex-conjugate regions to the 
upper and lower right of qinner, ^z,l,5 are dominant. In general, as the temperature T increases from to 
infinity for the Potts antiferromagnet, i.e. as a increases from to 1, the locus B moves in toward the origin 
in the q plane and degenerates to a point at the origin for a = 1. 

The boundary B q for the Potts ferromagnet is given for the illustrative finite-temperature value a — 2 in 
Fig. |ll]. Compared with analogous figures for the infinite-length limits of the L y — 2 strips of the square 
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Figure 8: Locus B q ({L}) for the Potts antiferromagnet on the m — > oo limit of the cyclic or Mobius L y = 2 strip of the 
honeycomb lattice for a = 0.1. For comparison, chromatic zeros are shown for the cyclic strip with m = 20 (i.e., n = 80). 
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Figure 9: Same as Fig. g, for a = 0.5. 
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Figure 10: Same as Fig. @, for a = 0.9. 
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Figure 11: Locus B q ({L}) for the Potts ferromagnet on the m — > oo limit of the cyclic or Mobius L v = 2 strip of the 
honeycomb lattice for a = 2. For comparison, chromatic zeros are shown for the cyclic strip with m = 20 (i.e., n = 80). 



17 



and triangular lattices (Figs. 13 in [|4j and flTql), one sees the common feature that B q passes through q = 0, 
does not cross the positive real q, but does have complex-conjugate lobes that extend into the Re(q) > 
half-plane. Similar to B q for this strip of the square lattice, here B q contains a finite line-segment on the 
negative q axis, whereas such a line segment is not present on B q for the triangular strip for this value of a. 

5.3 B in the Complex- Temperature Plane 

We next discuss slices of B in the complex-temperature plane, as a function of q. In Figs. [l^|l7| we show 
plots of B for q = 2, 3, 4, and 10. The locus B is noncompact in the a plane and, for q ^ 2, compact in the u 



plane. For q = 2 there can be noncommutativity of the type (1.9); the results that we show are for the f nq 
definition, i.e., we set q = 2 first before taking the limit n — > oo. For comparison with previous calculations 
of complex-temperature zeros in [p7[ p2[, we show some plots in the a plane as well as u plane. 
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Fi gure 12; Locus B U ({L\) for the m — > oo limit of the cyclic strip of the honeycomb lattice with q = 2. Partition function 
zeros are shown for m = 20, so that Z is a polynomial of degree e = 5m = 100 in v and hence, up to an overall factor, in u). 

In the u plane, four curves forming two branches of B cross at u = 0, reflecting the T — critical point of 
the g-state Potts ferromagnet. In general, the density of complex-temperature (i.e., Fisher zeros along 
the curves comprising B in the vicinity of a generic singular point u s behaves as |5p 

g-lu-usl 1 "** (5.3.1) 

where a s (a' s ) denotes the corresponding specific heat exponent for the approach to u s from within the CTE 
PM (FM) phase. Thus, for a continuous, second-order transition, with a s < 1, this density vanishes as 
one approaches the critical point u s along B, while, in contrast, if a s — 1, corresponding to a first-order 
transition, this density would remain nonzero as the boundary crosses u s , and there would be a discontinuity 
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Figure 13: Locus B U {{L}) for q = 3. Partition function zeros are shown for e = 100, as in Fig. 
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Figure 14: Locus £? a ({L}) for q = 3. Partition function zeros are shown for e = 100. 



19 



1 



0.5 



lm(u) o 



-0.5 



-1.5 -1 



-0.5 
Re(u) 



0.5 



Figure 15: Locus B U ({L}) for q = 4. Partition function zeros are shown for e = 100. 
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Fi gure 16: Locus B a ({L}) for q = 4. Partition function zeros are shown for e = 100. 
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Figure 17: Locus B U ({L}) for q = 10. Partition function zeros are shown for e = 100. 



in the internal energy, U . (We recall that this is analogous to the discontinuity in M{H) for T < T c as 
one changes H from positive through zero to negative values, and the connection to the nonzero density 
of Yang-Lee (complex-field) zeros, via the relation M(T) — 2ng(H — 0), where M(T) is the spontaneous 
magnetization |30|, |6l[.) Now the essential zero in the specific heat in (4.5) at T = 0, if expressed in terms 
of an algebraic specific exponent a, corresponds to a = — oo at tt = u s = 0. Substituting this into Q5.3.1 ), 
it follows that the density vanishes rapidly as one approaches the origin u = along the curves forming B u - 
This is evident in Figs. |lj and The angles at which these curves cross the origin are ±7r/4 and ±37r/4, 
as is clear in both of the figures shown, for q = 2 and 3. In particular, in the Ising case q = 2, because of 
the bipartite property of the cyclic strip of the honeycomb lattice, the antiferromagnetic and ferromagnetic 
versions of the model are equivalent, and so the antiferromagnet also has a T = critical point. This is 
manifested in the fact that B a passes through the origin of the a plane and B q passes through the point 
q = 2. Indeed, for this value q = 2, the locus B is the same in the a and u plane, i.e. it is invariant under 
the inversion map a — + a -1 = u. As was the case with the phase diagram for the infinite-length limit of 
the open strip, the positive u axis and its maximal complex-temperature analytic continuation forms the 
(CTE)PM phase. As discussed ear lier |L3|, fwj the fact that, for a = 0, the crossing of B q at q = 2 signals 
the property that the Ising antiferromagnet on the finite-width, infinite-length strips of the square lattice 
has a T = critical point is a useful feature of the periodic or reversed-orientation periodic longitudinal 
boundary conditions; in contrast, the loci B q at a = for strips with free longitudinal boundary conditions 
do not, in general, pass through q = 2. Furthermore, just as we showed by exact calculations in the case 
with the strips of the square lattice, in the present case of the honeycomb lattice, we anticipate that at the 
antiferromagnetic zero-temperature value a = 0, i.e., v = —1, the loci B q for wider cyclic/Mobius strips will 
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exhibit at least three crossings: (i) at a value q c that approaches the value for the infinite honeycomb lattice 

Qc(hc) = 3+ = 2.61803... (5.3.2) 



(the relevant root of eq. (5.3.3) for v — — 1, as discussed below), (ii) at q = 2, reflecting the T = critical 
point of the Ising antiferromagnet on these strips, and (iii) at q = 0. 

One of the reasons for interest in the complex-temperature phase diagrams for the g-state Potts model on 
infinite-length, finite-width strips of various lattices is that these provide exact results that exhibit features 
which can give insight into the properties of the complex-temperature phase diagrams of the corresponding 
Potts model on the respective 2D lattices jlj, Of course for the q = 2 Ising case, since the model is 
solved, this complex-temperature phase diagram is known J5(J. However, for other values of q, these insights 
can help one in interpreting the information that one has from calculations of complex-temperature zeros 
on finite sections of the honeycomb lattice |5?j, ^32). In order to appreciate some of the similar features in 
the complex-temperature phase diagrams, we recall that B a for the Ising model on the honeycomb lattice 
consists of an arc of the unit circle a = e ±l8 with 7r/3 < 9 < tt, together with a closed, bean-shaped 
curve that intersects the real axis at a = 2 — -\/3, the paramagnetic-antiferromagnetic critical point, and its 
inverse, a = 2 + \/3, the paramagnetic-ferromagnetic critical point [Q. The arc of the unit circle and the 
closed curve intersect at a = ±i. There are, of course, obvious qualitative differences between the locus B a 
for the finite-width, infinite-length strip and the full 2D honeycomb lattice, reflecting the fact that the 2D 
Ising ferromagnetic (and equivalently, the antiferromagnet) has a finite-temperature critical point, while this 
occurs only at zero-temperature on the infinite-length, finite-width strip, which is effectively a quasi-one- 
dimensional system. However, for both the strip and the full honeycomb lattice, B a includes intersection 
points where two branches of curves cross at a = ±i. The locus B a for the infinite- length L y = 2 strip also 
exhibits complex-conjugate arc endpoints in the Re(a) > half-plane, as is true of the exactly solved Ising 
model on the full 2D honeycomb lattice; in addition, this locus for the strip exhibits complex-conjugate arc 
endpoints in the Re(a) < half-plane. This illustrates how, although physical temperature properties of the 
Potts model on an finite- width, infinite strip of a given lattice and on the full 2D lattice are clearly different, 
complex-temperature properties can have some similarities. 

Having discussed the situation with the q = 2 case, where the 2D model is exactly solved, let us proceed 
to cases where the 2D Potts model has not been solved. In general, for the honeycomb lattice, from duality 
and a star-triangle relation, one can derive an algebraic equation for the PM-FM and PM-AFM critical 

points §, m n 

v 3 -3qv-q 2 = . (5.3.3) 

(where v = a — 1). For example, for q = 3, the equation yields three solutions (tabulated, e.g., as eqs. 
(3.3)-(3.5) in §7|), namely (i) the PM-FM critical point at a PM -FM, q =3 = 1 + 2\/3 cos(vr/18) = 4.41147..., 
and two complex-temperature singular points, at a = 1 — V3 cos(7r/18) + 3sin(7r/18) = —0.1847925... and 
a = 1 - a/3cos(7t/18) - 3sin(7r/18) = -1.22668.... There is no physical PM-AFM critical point solution, 
and the absence of such a phase transition was confirmed by explicit Monte Carlo study of the model in 
]54j (and further confirmed in |55|). More generally, considering this solution as a function of q generalized 
from positive integer to real values, one sees that it decreases from the q — 2 Ising value a = \[2 — 1 and 



passes through a = at the value q c (hc) given in eq. (5.3.2) above (see the middle curve in Fig. 4 of 
f37f), corresponding to a zero-temperature critical point for the Potts antiferromagnet at this value of q. For 
q c {hc) < q < 4, this root decreases from to —1, where it coalesces with another solution (the lower curve in 
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Fig. 4 of j57|) to form a double root, and for q > 4, eq. ( 5.3.3j ) has only a single real root, namely the PM-FM 
critical point. In [5(1 , as a rigorous result of duality, it was shown that as a consequence of the (first-order) 
phase transition of the q = 3 Potts antiferromagnet on the triangular lattice at apM-AFM.t — 0.203, it 
follows that for q = 3 Potts model on the honeycomb lattice, the left-most crossing of the locus B a is at the 
complex-temperature point given by the duality mapping v — ► q/v, i.e., 

a — 1 

i.e., a ~ —2.76. The calculations of complex-temperature zeros in []57| (given in Figs. 5-7 of that paper) were 
consistent with these crossings of B a at a ~ —2.76, —1.23, and —0.185 and also suggesteda fourth possible 
crossing, at a ~ —0.65, for a total of four crossings on the negative real a axis. 

For the q = 3 Potts model on the infinite-length strip, as shown in Fig. we find that B a crosses 
the negative real a axis at two points, a = —1/2 and a = —2. In order to relate this complex-temperature 
phase diagram to one for the q = 3 Potts model on the 2D honeycomb lattice, one can imagine performing 
a transformation in which one retracts the branches of the curves forming B a from the zero-temperature 
ferromagnet critical point at u — 0; in the right-hand half-plane, one would then reconnect the two branches 
of the curve so as to cross the positive real a axis at the value of a,pM-FM,q=3 given above, while in the 
left-hand half plane, one would reconnect the branches to cross at a point on the negative real a axis. One 
could also pull apart the two curves that have a tacnode multiple point at a = —1/2 to obtain a total of four 
crossings on the negative real axis, as appears to be the case for the model on the full 2D honeycomb lattice. 
The exactly determined locus B a also exhibits a complex-conjugate pair of arc endpoints in the Re(a) > 
half-plane, qualitatively similar to the arc endpoints that appear to be formed by the complex-temperature 
zeros for 2D honeycomb lattice. Thus, although the physical thermodynamic properties are clearly different, 
certain complex-temperature features of the phase diagram for the q — 3 Potts model on the infinite-length, 
width L y = 2 strip of the honeycomb lattice show interesting qualitative similarities to those inferred for the 
complex-temperature phase diagram of the model on the full honeycomb lattice. 

We next consider the case q — 4. For this value, eq. ( 5.3.3 ) reduces, in terms of the variable a, to 
(a — 5) (a + l) 2 = 0, yielding the physical PM-FM critical point apM-FM,q=4 — 5 and a double root at the 
complex-temperature point a = u = — 1. As in the q — 3 case, there is no phase transition at finite or zero 
temperature for the q = 4 Potts antiferromagnet. From the duality relation given in [[56| together with the 
fact that the q = 4 Potts antiferromagnet has a T = critical point on the triangular lattice, it follows that 
the left-most crossing of B a on the real axis for the q = 4 Potts model on the honeycomb lattice is at the 
value of a given by ( 5.3.4 ), viz., a = —3. The calculations of complex-temperature zeros in |62|| (see Figs. 
2-4 of that paper) exhibited these three crossings of the complex-temperature phase boundary B a at a = 5, 
a = — 1 and a = — 3 and also suggested crossings at approximately a — — 2 and a — —1/2. The locus B a 
also exhibited a pair of arc endpoints (prongs) in the Re{a) > half-plane and structures that resembled 
a complex-conjugate pair of prongs or bulbs in the Re(a) < half-plane, all within the interior of the PM 
phase. 

What insights can we draw for this q = 4 case from our exact results for the infinite-length, L y = 2 
strip? The locus B is shown in the u and a planes in Figs, [l^ and [l6]. In the a plane, we again perform a 
transformation to remove the intrinsically ID features of the complex-temperature phase boundary, i.e. their 
passage through u = for the comparison with the 2D results. We thus retract the branches of the curves 
passing through this point; in the right-hand half-plane, we reconnect the two branches of the curve so as 
to cross the positive real a axis at the value of apM-FM,q=i given above, while in the left-hand half plane, 
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we reconnect the branches to cross at a point on the negative real a axis, which then forms the left-most 
crossing. We could also pull apart the curves that meet and cross at a = — 1, thereby obtaining four crossings 
of B a on the negative real a axis. The two complex-conjugate pairs of arc endpoints on B a are in nice 1-1 
correspondence with those observed for the model on the 2D honeycomb lattice in |6^|. In addition, B a 
exhibits two small bulb-like structures around a = — 0.65±0.53i; these bear an interesting correspondence to 
the analogous structures found in ]6^| . This discussion shows how some qualitative aspects of the complex- 
temperature phase diagram for the Potts model on the 2D honeycomb lattice can correspond to those found 
via exact solutions for infinite-length, finite-width strips. One of the most remarkable findings is that the 
crossing point of B a at a = u = — 1 for the infinite-length strip is not just qualitatively but exactly the same 



as one finds, rigorously, from the solution of (5.3.3), for B a on the infinite honeycomb lattice. This is also 
true for the infinite triangular lattice fl62[ , since it is dual to the honeycomb lattice, and under the duality 
mapping ( |5.3.4 ), a = — 1 maps to itself if q = 4. As was discussed in ^35|, the point a = —1 is a solution of 



the equation for the critical manifold, v 2 = q, at q = 4 and hence B a passes through a = — 1 for the square 
lattice. Calculations of complex-temperature zeros for the q = 4 Potts model on large finite sections of the 
triangular lattice in (6^] and the square lattice (see |64|, ^33[ |65| and references therein) are both consistent 
with these exact results. Moreover, our exact determinations of B a for the ID Potts model (i.e., n — > oo 
limit of the circuit graph) ]5ll |TJ] and infinite-length, L y — 2 cyclic/Mobius and open strips of the square 
and triangular lattices (IJ, [l6| showed that in all cases, if q = 4, then B a passes through a = — 1. It is of 
interest to see what the values 5 = 4, a = v + 1 = — 1 correspond to in terms of the variables x and y, defined 



in eqs. ( 13.2.3 ) and ( 13.2.4 ), that enter in the Tutte polynomial equivalent to the Potts model partition 



function (see eqs. (13.2.1) and (13.2. 6|) ) . We find that (q,v) = (4,-2) corresponds to the symmetric point 



(x, y) = (—1, —1). More generally, one can investigate the continuous accumulation set of zeros of the Tutte 
polynomial in the complex x plane for fixed y, B x and in the complex y plane for fixed x, B y . Since this 
information is equivalent to our study of the loci B q for fixed u and B u for fixed q, we shall not include the 
results here. 

An illustration of B u for larger values of q is given for the case q = 10 in Fig. [l7| This locus includes 
two separate figure-eight type curves, the one to the left also involving outgoing arcs, and, to the far left, a 
self-conjugate arc. As compared with our result for the infinite-length L y — 2 strip of the square lattice in 
Fig. 17 of |Q, one sees that the middle figure-eight curve plus outlying arcs that are present on B u for the 
honeycomb strip are absent in the locus for the strip of the square lattice, while the other two components 
(the self-conjugate arc on the left and the figure-eight curve on the right) are qualitatively similar. For the 
infinite-length honeycomb lattice strip with open boundary conditions, we note that B u consists of the same 
self-conjugate arc on the left as in Fig. |l7], together with two complex-conjugate pairs of arcs. 

6 L y = 3 Cyclic and Mobius Strips of the Honeycomb Lattice 

Illustrations of strips of the honeycomb lattice with width L y — 3 and (FBC y , PBC X ) — cyclic and 
(FBCy, TPBC X ) = Mobius boundary conditions, displayed as strips of a brick lattice, are shown in Fig. [l8| 
(a,b). 
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Figure 18: Illustrative strip graphs of the honeycomb lattice, displayed as a brick lattice, with width L y = 3: (a) cyclic, 
length m = 3 and (b) Mobius, length m = 2. We also show (c) an L y = 3 cyclic crossing-subgraph strip of the honeycomb 
lattice with length £ = 5. As discussed in the text, the cyclic crossing-subgraph strips of a given width arc equivalent to the 
cyclic and Mobius strips of the same width for even and odd £, respectively. This equivalence is illustrated here; graph (c) is 
equivalent to (b). 



We calculate the chromatic polynomials by iterated use of the deletion-contraction theorem and by 
coloring matrix methods. Here and below, in our results for chromatic polynomials, we shall suppress the 
subscript P in \p,G,j an d cp.G.j where it is obvious. For the L y = 3 cyclic and Mobius strips of the honeycomb 
lattice, denoted hc3c — he, 3 x m, FBC y , PBC X and hc3mb — he, 3 x m, FBC y , TPBC X , we obtain the result 
N h c3c.x = N hc3mbA = 14 and 



where, for j = 1, 2, 



and 



Note that 



14 



P(hc, 3 x m, FBC y ,PBC x ,q) = ^ c hc3C}j {\ hc3J ) m 

3=1 

14 

P{hc,3 x m,FBC y ,TPBC x ,q) = ^ c hc3mbiJ (Xhcs.jY 

3=1 

A/ic3,l = 1 

Ah c3 ,2 = {q- 1) 2 • 

^hc3j = ^hc2, 3 for j = 1, 2 . 



(6.1) 

(6.2) 

(6.3) 
(6.4) 
(6.5) 
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The \hc3,j for j — 3, 4, 5 are the roots of the cubic equation 

£ 3 - (2q 2 -8q + 11)£ 2 + (g 4 - 8<? 3 + 26q 2 - 38q + 23)£ - (<? - l) 2 = (6.6) 
and the \hc3,j f° r 6 < J < 11 are the roots of the sixth degree equation 
C 6 - (3g 4 - I8q 3 + A6q 2 - 60q + 37)£ 5 

+ {3q 8 - 36q 7 + 199g 6 - 662q 5 + 1456g 4 - 2172g 3 + 2144g 2 - 1278g + 358)£ 4 
-{q 12 - 18q n + lMq 10 - 828q 9 + 3117q s - 8646g 7 + 18070? 6 - 28564g 5 + 33782g 4 
-29110.? 3 + 17370 g 2 - 6466g + 1142)£ 3 

(q - lf{q 12 - 20q u + I86q w - 1060g 9 + 4115g 8 - 11454g 7 + 23448g 6 

-35642 9 5 + 40094 g 4 - 32780 9 3 + 18680g 2 - 6756g + 1197)£ 2 

-(q - l) 4 {2q 8 - 26/ + 150g 6 - 498g 5 + 1037g 4 - 1388 9 3 + 1176g 2 - 590q + 141)£ 

+ (g-2) 2 ( g -l) 8 = 0. (6.7) 
Finally, \hc3,j for 12 < j < 14 are roots of the cubic equation 

£ 3 - (q 6 - 8q 5 + 28q 4 - 56q 3 + 71q 2 - 58q + 26)£ 2 

+ (q - l) 2 (q 6 - 10q 5 + 43g 4 - 102g 3 + 144g 2 - 120g + 49)£ 

-(q-2) 2 (q-l) i = 0. (6.8) 

The corresponding coefficients for the cyclic strip are 

ChcScA = c (3) (6.9) 

c h c3c.j = c (2) for 2 < j < 5 (6.10) 

c hc 3c,j = c (1) for 6 < j < 11 (6.11) 

ChcSc,j = c<°> = 1 for 12 < j < 14 . (6.12) 

For the Mobius strip, as will be discussed further below in section [t], the coefficients are not, in general, 
simply ±1 times Chebyshev polynomials of the second kind, c^ d \ but instead are products of the with 
certain terms Xcghcs.j's occurring in the chromatic polynomial for related families of strip graphs called 



crossing-subgraph strips, defined via eqs. (|7.2j )- (|7.8| ). We have 

c hc3m b,i = -c (2) = -(q 2 -3 g +l) (6.13) 

Chc3mb,2 = -C (0) A cff ^ c3 ,3 = ~(q - 1) (6-14) 
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Chc3mbj ■— -c (0) A C gh c3 j + 2 for j = 3, 4, 5 (6.15) 

Chc3mb,j = c (1) A C5/lc3j+ 5 for 6 < j < 11 (6.16) 

Chc3mb,j = c (0) A c ^ C 3j +5 for 12 < j < 14 . (6-17) 
By general coloring matrix methods, we have 

Ly-l 

Cp,G = q(q(q-V)) 2 for G = hc,L y x m,FBC y ,PBC x , L y odd (6.18) 

Cp,g = (q{q - 1)) 2 for G = he, L y x to, FBC y ,PBC x , L y even (6.19) 
Our results for the Mobius strips lead to the inference that 

Cp, G = for G = he, L y x to, FBC y , TPBC X . (6.20) 

The first special case of this, for L y — 2, was obtained in p8[ . 




Fi gure 19; Chromatic zeros for the Ly = 3 cyclic strip of the honeycomb lattice of length m = 20, i.e., n = 120 vertices. 

As a special case of the arguments given in Q, since the Agj 's are the same for the cyclic and Mobius 
strips, and since in the limit m — * oo the continuous singular locus B is defined as the solution to the 
degeneracy of leading terms Xgj, it follows that this locus is the same for the cyclic and Mobius strips. In 
Fig. [l9| we show the chromatic zeros for a cyclic strip of the honeycomb lattice with width L y = 3 and 
length, to = 20, i.e., n = 120 vertices. This length is sufficiently great that the chromatic zeros for the cyclic 
and Mobius strips are generally similar and lie close to the asymptotic singular locus B. We find the exact 
result 

q c = 2 for he, 3 x oo, FBC y , (T)PBC X . (6.21) 



27 



This is the same value that was obtained in |38) for the (m — > oo limit of the) L y — 2 cyclic or Mobius strip 
of the honeycomb lattice (below we shall show that the L y = 4 strip with cylindrical boundary conditions 



yields a value of q c closer to the value (5.3.2) for the infinite 2D honeycomb lattice). The locus B for L y = 3 
has support for Re(q) < as well as Re(q) > 0, just as was true for L y = 2 |38[ . 

This L y = 3 locus separates the q plane into several regions. The outermost one, region R\, extends to 
infinite |g| and includes the intervals q > 2 and q < on the real q axis. Region i?2 includes the real interval 
< q < 2. Two complex-conjugate regions R3,R%, centered at approximately q = 2.25 ± 0.55i, occur to 
the upper and lower right of q c , as is evident in Fig. |l^. Two further complex-conjugate regions, R^,R^, 
centered at approximately q = 1.9 ± 0.5i, occur to the upper and lower left of q c . One can observe that, 
just as was the case for L y = 2, the density of zeros is small on the boundary separating region R4 from 
i?2 and its complex conjugate boundary separating R\ from i?2- As previous exact calculations have shown 
p6| , pTU| , fill , region diagrams for infinite- length, finite- width strips of lattices can also include extremely small 
regions. We have found evidence for a complex-conjugate pair of these tiny regions to the upper and lower 
right of q c . We have not made an exhaustive search for other tiny regions. 

In region Ri, the dominant \hc3,j is the root of the cubic equation (6.8) with the largest magnitude, 
which we denote as Xhc3,i2- Hence, 

W = (X hc3A 2) 1/6 , qeRx. (6.22) 

The fact that this is the same as W for the (FBCj,,FBC x ) case is a general result | fL3| . In region R2, 
and also in regions Rs,R^, the dominant terms are the roots of eq. (6.7) that are largest in magnitude in 
these respective regions. We shall refer to these respective dominant terms generically as \hc3fi, so that 

\W\ = |A, ic3 , 6 | 1/6 , qeR 2 ,R 3 ,R* 3 . (6.23) 



In regions i?4,i?|, the dominant term is the root of eq. ( J6 . 6[ ) with maximal magnitude, which we denote as 
Aftc3,3, so that 

\W\ = |A ftc3 , 3 | 1/6 , qeR 4 ,Rt. (6.24) 

Note that W = 1 at q c = 2. A summary of some features of these strip graphs is given in Table [lj Regarding 
the comparison with previous results, those for the L y = 2 cyclic/Mobius strips are from Ref. |}8) and those 
for the L y = 3 open strip are from Ref. p5| . 

7 L y = 3 Cyclic Crossing- Subgraph Strips of the Honeycomb Lat- 
tice 

An illustration of the cyclic L y = 3 crossing-subgraph strip of the honeycomb lattice, denoted cghc3, of length 
t = 5, is shown in Fig. [l8| (c). For this type of strip, with arbitrarily great £, we calculate N cg h C 3. C yc. = 19 
and 

19 

P{cghc, 3xf, FBC y , PBC X , q) = ^ 

CcghcS.j 

where for 1 < j < 4, 

X C ghc3,l = — A cff /ic3.2 = 1 (7-2) 
X C ghc3,3 = —X cg hc3A = 1 — 1 • (7-3) 
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Table 1: Properties of P, W , and B for strip graphs G s of the honeycomb (he) lattice. The properties apply for a given strip 
of size L y X m; some apply for arbitrary m, such as N\, while others apply for the infinite-length limit, such as the properties 
of the locus B. For the boundary conditions in the y and x directions (BC y , BC X ), F, P, and T denote free, periodic, and 
orientation-reversed (twisted) periodic, and the notation (T)P means that the results apply for either periodic or orientation- 
reversed periodic. The column denoted eqs. describes the numbers and degrees of the algebraic equations giving the Ag s j ; for 
example, {2(1), 2(3), 1(6)} indicates that there are 2 linear equations, 2 cubic equations and one sixth degree equation. The 
column denoted BCR lists the points at which B crosses the real q axis and the value of q c for the given family. The notation 
"none" in this column indicates that B does not cross the real q axis. The notation "int;qi, q c " refers to cases where B contains 
a real interval, there is a crossing at qi, and the right-hand endpoint of the interval is q c . Column labelled "SN" refers to 
whether B has support for negative Re(q), indicated as yes (y) or no (n). 



Ly 


BCy 


BC X 


N x 


eqs. 


BCR 


SN 


2 


F 


F 


1 


{1(1)} 


none 


n 


3 


F 


F 


3 


{1(3)} 


2 


n 


4 


F 


F 


5 


{1(5)} 


int;2.093,2.099 


n 


2 


F 


(T)P 


4 


{4(1)} 


0, 2 


y 


3 


F 


(T)P 


14 


{2(1),2(3),1(6)} 


0, 2 


y 


4 


P 


F 


4 


{1(4)} 


int;2. 222,2. 250 


y 



The X cg hc3.j for j = 5, 6, 7 are the roots of the cubic equation 

e+Z 2 -(q 2 ~M + 5)t~(q-l) = ■ (7.4) 
The X cg hc3,j for j — 8, 9, 10 are the roots of the cubic equation 

e-e-(q 2 ~M + ^+(q-l) = . (7.5) 



Since eq. ( |7.5| ) is related to eq. (7^4) by the reversal of the signs of the coefficients of the £ 2 and £° terms, 



it follows that the respective roots of these equations are opposite in sign: 

Kghcz.j = -X cg hc3,]+3 for j = 5, 6, 7 . (7.6) 
The X C ghc3,j for 11 < j < 16 are the roots of the sixth degree equation 

C 6 + (q 2 -M + 5)£ 5 - (g 4 - 5q 3 + 10q 2 - lOq + 6)£ 4 

-{q 6 - 10q 5 + 43q 4 - 102q 3 + U2q 2 - 110g + 38)£ 3 

+ {q - l)(q 6 - 10q 5 + 43g 4 - lOOg 3 + 132g 2 - 94q + 29)£ 2 

+(q - 1) V -Aq + 5)£ -(q- 2)(q - l) 4 - . (7.7) 
Finally, the X cg hc3,j for j — 17, 18, 19 are the roots of the cubic equation 

f - (q - 2){q 2 -2q + 3)£ 2 + (q - l)(q 3 - 5q 2 + 8q - 5)£ - (q - 2) (q - l) 2 = . (7.8) 
The corresponding coefficients are 

Ccghc3,l - i(c (3) - c (2) ) = ~{q- 2)(q 2 -4q + l) (7.9) 
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c cghc3 . 2 = l(c {3) + c (2) ) = \q{q - l)(g - 3) (7.10) 

c cghc ^ = \(c {2) - C (°)) = ig(g-3) forj = 3,5,6,7 (7.11) 

c cg h c3j = ^(c (2) +c(°)) = i( ( 7-l)( ( 7-2) for j = 4, 8, 9, 10 (7.12) 
c C ghc3,j = q - 1 for 11 < j < 16 (7.13) 
Ccghcsj = 1 for 17 < j < 19 ■ (7.14) 

In jl8| we have discussed how cyclic crossing subgraph strips of the square lattice with free transverse 
boundary conditions reduce, for even and odd length L Xl to cyclic and Mobius strips, respectively. Similarly, 
here, for even and odd length £, the cyclic crossing subgraph strip of the honeycomb lattice reduces to the 
cyclic and Mobius strips of this lattice, so that the chromatic polynomial for this cyclic crossing subgraph 
strip reduces to the respective chromatic polynomial of the cyclic and Mobius strip. With our labelling 
conventions for the length of the cyclic crossing subgraph strips as compared with the cyclic or Mobius strips 
of the honeycomb lattice, the cyclic crossing subgraph of even length £ = 2m (odd length I = 2m + 1) is 
identical to the cyclic (Mobius) strip of this lattice of length to, respectively. There is a resultant reduction 
in the number of Ac.j's, from 19 to 14, as a result of the fact that there are five pairs of terms X cg hc3,j 
that differ in sign and have different coefficients, as given above. For the even-length case we thus have the 
identifications 

(K g hc3j) 2m = (>~hc3j) m (7-15) 



while for the odd-length case, 



{X C ghc3,j) m — ^cghc3,j(^hc3,j) m (7-16) 



so that, for both the even- and odd-£ cases for the crossing-subgraph strips, i.e., for both the cyclic and 
Mobius strips, 

A^ c3 ,i = (A Cff ^ c3j ) 2 for j = 1,2 (7.17) 

Xhc3,2 = (A c ^ c3j ) 2 for. ] -3, 4 (7.18) 

Xhc3,3 = (Kghc3,j) 2 for. ] -5, 8 (7.19) 

Afccs.4 = (A c ^ c3j ) 2 for. j = 6, 9 (7.20) 

Xhc3.5 = (Xcghc3.j) 2 for j = 7, 10 (7.21) 

Ahcsj = (A cg/lc3j+ 5) 2 for 6 < j < 14 . (7.22) 



In the case of odd length, we absorb the left-over factor of X cg hc3.j in eq. (7.16) in the definition of the 
coefficients. The relations (7.17)-(7.22) expressing A^ c3j 's for the cyclic/Mobius strips as squares of A cg h C 3 ; fc's 
are reflected in obvious relations between the algebraic equations for these terms. 

In the limit £ — * oo, one derives the singular locus B for the crossing-subgraph strip. As discussed in 
p8| , since one can take this limit on even values of £ or, alternatively, on odd values of £, it follows that the 
singular loci B for the infinite- length limits of the cyclic and Mobius strips of the honeycomb lattice are the 
same. Here this already follows from the property that the Agj 's are the same for the cyclic and Mobius 
strips. 
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8 Structural Properties of Chromatic Polynomials for Cyclic Strips 
of the Honeycomb Strip 

Let us define np(hc, L y ,d) as the number of terms Xp,hc,L y ,j, each of which has, as its coefficient, the 
Chebyshev polynomial cS d \ in the chromatic polynomial P(hc,L y ,q) of the cyclic strip of the honeycomb 
lattice of width L y , denoted hc,L y . Clearly, the total number, Np^ c ^ \, of different terms Xp t hc,L y .j in 



(1.13) for this family of strip graphs is given by 

Ly 

A ;<./„. /...a = ^2 n p{hc, L y ,d) . (8.1) 



d=0 



Using the general formulas ( 6.1 8j ) and ( 6.19| ) together with the coloring matrix methods that we have discussed 



in detail in p7[ , we can calculate the numbers np(hc, L y , d) and Np,h c ,L y ,\- We find 

n P (hc,L y ,d) = Q for d > L y , (8.2) 

n P (hc, L y , L y ) = 1 (8.3) 

and the values np(hc, 2,0) = 1, np{hc, 2,1) = 2. All other numbers np(hc, L y ,d) are then determined by 
the following recursion relations. For even L y > 4 

rip (he, L y , d) — np(hc, L y — 1, d — 1) + ripihe, L y — 1, d) + np(hc, L y — 1, d + 1) for d > 1 (8-4) 

while for d = 0, 



np(hc, L y , 0) = np(hc, L y — 1,1) . (8.5) 



For odd L y > 3 



np(hc, L y ,d) = np(hc, L y — 1, d — 1) + 2np(hc, L y — l,d) + np(hc, L y — 1, d + 1) for d > 1 (8.6) 
while for d = 0, 

np(hc, L y , 0) = np(hc, L y — 1,0) + np(hc, L y — 1,1) . (8-7) 



We observe that the recursion relations (p.4|) and (8.5) are the same as we found in [|47j for the numbers 



np(L y ,d) for cyclic strips of the square and triangular lattices, while eqs. (8.6) and (3.7) are analogous to 
the recursion relations that we found for the numbers nz(L y , d) for cyclic strips of the square and triangular 
lattices (in JitJ these applied for both even and odd L y ). We also obtain the relations 

N PM ,L y ,X = 4JV> 1 h c ,L B _i 1 A - 2n P (hc, L y - 1, 0) for odd L y > 3 (8.8) 

and 

Np t hc,L y ,x = SNp,hc,L v -i,x - 2n P (hc, L y - 1, 0) for even L y > 4 (8.9) 

whence 

N P>hc ,L y ,\ = l2N PMtLy ^ x - 8n P (hc, L y - 2, 0) - 2n P (hc, L y - 2, 1) for L y > 4 . (8.10) 

Results for L y up to 10 are given in Table |^. These may be compared with the analogous calculations 
that we carried out for cyclic strips of the square and triangular lattices in fl47| . For example, the values 
Np : hc,L y ,\ = 4, 14 and 36 for L y = 2, 3 and 4 may be contrasted with the numbers Np t a,L v ,\ = 4, 10 and 26 
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Table 2: Table of numbers np(hc, L y , d) and their sums, Np t hc,L y ,X f° r cyclic strips of the honeycomb lattice. Blank entries 
are zero. 



L y i d — > 
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8 


352 
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3114 
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321 


81 


13 


1 




14288 


10 


3114 


8025 


9840 


8625 


5847 


3123 


1305 


415 


95 


14 


1 


40404 



for G — sq, tri with the same set of L y values. It is straightforward to extend these to higher values of L y , 
and related results can also be given for Mobius strips of the honeycomb lattice. 

From eq. ( |B.1C ), it follows that the leading asymptotic behavior of Np^ c ,L y ,x (ignoring power-law pref- 
actors) is 

N P , hc ,L vl x ~ (12) L «/ 2 = (3.4641. ..) L y as L y -> oo . (8.11) 

This may be contrasted with the leading asymptotic behavior of the corresponding numbers for the cyclic/Mobius 
strips of the square and triangular lattices |^7j] 

Np,G,L y ,\ ~ 3 L " for G = sq 7 tri as L y -> oo . (8.12) 



Regarding the full Potts model partition function, we note that a generalization of the relation (5.1.15) 
to Cz = q Ly for the cyclic honeycomb strip with width L y yields the result N Zi hc.L y ,\ — ( 2 l")j similarly to 
the case with the corresponding strips of the square and triangular lattices. 



9 Cylindrical L y = 4 Strip of the Honeycomb Lattice 

The first width for which a strip of the honeycomb lattice with cylindrical = (PBC y , FBC X ) boundary 
conditions can naturally be defined, satisfying A max = 3, is the width L y = 4. A picture of such a lattice 
is shown in Fig. ^ for length m = 3 bricks. In the labelling convention of |3f| [?7|, this would be denoted 
as m = 2; the present convention yields somewhat simpler expressions for the coefficient functions in the 
numerator of the generating function. 



32 





> < 


1 < 


• 1 


> < 


| ( 


































> H 


» r< 


1 r« 


1 H 


> — 


» H» 



Figure 20: Illustrative strip graph of the honeycomb lattice with width L y = 4, length m = 3 bricks and (PBC y , FBC X ) = 
cylindrical boundary conditions. 

In general, for the strip of the honeycomb lattice with cylindrical = (PBC y , FBC X ) boundary conditions 
there are n = 2L y (m + 1) vertices and e = L y (3m + 2) edges. As before, it is convenient to present the 
results in terms of a generating function, denoted T(G s ,q,x), where G s refers to the type of strip graph. 
The chromatic polynomial P((G s ) m ,q) is determined as the coefficient in a Taylor series expansion of this 
generating function in an auxiliary variable z about z = 0: 



T(G s ,q,z)= ^P((G s ) m , g ) 2 



m=0 



The generating function T(G s ,q, z) is a rational function of the form 

jV(G s ,q,z) 



T(G s ,q,z) = 



V(G s ,q,z) 



with 



and 



M(G s ,q,z) = J2 A G e ,Mz J 

3=0 



V(G s ,q,z) = l + J2l>G a A<l> j 

3 = 1 



(9.1) 



(9.2) 



(9.3) 



(9.4) 



where the Aq 8 j and bc s j are polynomials in g, and djsf = deg z {N), dx> = deg z (V), In factorized form 



V(G s ,q,z) = Y[(l-\ GsJ (q)z) . 

3 = 1 



(9.5) 



We find that deg z (j\f) = 3 and deg z (T>) — N\ = 4. The coefficient functions in the numerator and 
denominator of the generating function are 

A = q(q - l)(q 6 - 7q 5 + 2lq A - 35g 3 + 35q 2 - 2lq + 7) (9.6) 

A 1 = -q(q- l){q 10 - I5q 9 + 103q 8 - 417q 7 + 1074g 6 - 1782g 5 + 1850g 4 

-1125<7 3 + 422g 2 - 246(7+ 161) (9.7) 
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A 2 = g (g-l)(V -6V + 460 9 8 - 1956 9 7 + 5461g 6 - 10497<7 5 

+ 14128g 4 - 13212q 3 + 8281g 2 - 3203<? + 614) (9.8) 

As = -4q{q-l) B (q-2) 2 (9.9) 

b x = -q 8 + I2q 7 - 66q 6 + 220g 5 - 496g 4 + 796<? 3 - 922g 2 + 734g - 314 (9.10) 

b 2 = 9 12 -209 n + 188(? 10 - 1092g 9 + 4344q 8 - 12428g 7 + 26192g 6 

-41022g 5 + 47597<? 4 - 40318g 3 + 24247<? 2 - 9782g + 2169 (9.11) 

6 3 = -4? 12 + SOq 11 - 736g 10 + 4124g 9 - 15705<j 8 + 42922g 7 - 86535g 6 

+129964g 5 - 144575g 4 + 116374g 3 - 64525g 2 + 22280g - 3680 (9.12) 

b 4 = A[{q-l){q-2)f . (9.13) 
The W function is given as W = (Xhcipfj) 1 ^ where \hdpf,j is the maximum-magnitude root of the equation 

£ 4 + &i£ 3 + 6 2 £ 2 + ^ + &4 = (9.14) 



with the bj coefficients given in eqs. (|9~10|)-(|9~ll) 



A plot of the chromatic zeros for a long finite strip of this type is shown in Fig. |2l]. For this great a length, 
these zeros give a reasonably accurate indication of the position of the asymptotic locus B. From our exact 
analytic results, we find that the locus B includes arcs and a short interval on the real axis 2.194004 < q < q c , 
where q c = 2.249840. It is interesting that this value of q c is just 14 % below the value ( 5.3.2] ) for the infinite 



2D honeycomb lattice. We also note that the locus B includes support for Re(q) < 0, as is evident from Fig. 

10 Free L y = 4 Strip of the Honeycomb Lattice 

In general for a strip of the honeycomb lattice with free boundary conditions, width L y , and m bricks, there 
are n = 2(L y m + L y — 1) vertices and e = (3L y — l)m + 2L y — 3 edges. An illustration of the strip of the 
honeycomb lattice with width L y = 4, length m = 3, and free = (FBC y , FBC X ) boundary conditions is 
shown in Fig. (in the labelling convention of |33, pTi], this would be denoted as length m = 2). 
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Figure 21: Chromatic zeros for the width L y = 4 strip of the honeycomb lattice with cylindrical boundary conditions, 
(PBC y ,FBC x ), for length m = 14 (i.e., n = 120). 



Figure 22: Illustrative strip graph of the honeycomb lattice with width L y = 4, length m = 3 bricks and open boundary 
conditions. 
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We find that deg x (N) = 4 and deg x (T>) = N\ = 5. The coefficient functions in the numberator and 
denominator of the generating function are 

A Q = q(q-lf (fO.f) 
Ai = -q(q - 1)(2/ - 22q 7 + 109/ - 314/ + 566/ - 636/ + 41 7q 2 - 133g + 10) (10.2) 

A 2 = q(q- l/O? 11 - 17/° + 133/ - 628/ + 1975/ - 4306/ + 6570/ 

-6933/ + 4921/ - 2273/ + 683<? - 133) (10.3) 

A 3 = -q(q- l) 2 (2q n -34/° + 268/ - 1289/ + 4182/ - 9547 / 

+15548/ - 17985/ + 14487/ - 7811/ + 2599g - 424) (10.4) 
A 4 = 2q{q - l) 5 (q - 2){2q 5 - 16/ + 49/ - 7l/ + A9q - 12) (10.5) 

hi = -{q 2 -Aq + 5)(/ - 7q 5 + 22q 4 - 42c? 3 + 54/ - 47q + 27) (10.6) 
b 2 = 2q 12 - 36c/ 1 + 305/° - 1609/ + 5894/ - 15827c/ 7 + 32046/ 

-49513/ + 58255c/ - 51198t7 3 + 32069/ - 129 1 5c? + 2555 (10.7) 

b 3 = -/ 6 + 25/ 5 - 297/ 4 + 2224/ 3 - 11745/ 2 + 46379/ 1 

-141680/° + 341694c7 9 - 658019c7 8 + 1016383c7 7 - 125 6 9 55g 6 

+1233979/ - 945308c7 4 + 54 8 4 23c7 3 - 228538c7 2 + 6 1 589g - 8161 (10.8) 

6 4 = (c7-l) 2 (2/ 4 - 44/ 3 + 456/ 2 - 295l/ 1 + 13318/° -44313/ 

+112043c7 8 - 218692/ + 331398c/ - 388712/ + 348419/ 

-232363/ + 109612/ - 33033(7 + 4876) (10.9) 

b 5 = -2(q - l) 6 (c7 - 2) 2 (2/ - 20/ + 83/ - 185/ + 239/ - I7bq + 60) . (10.10) 
The W function is given as W = (^hcif f.j) 1 ^ 8 , where \hdff,j is the maximum- magnitude root of the equation 

e + &i£ 4 + b 2 e + b 3 e + & 4 $ + & 5 = o ao.n) 



with the bj coefficients given in eqs. ( 10.6 )-( 10.1C ). 

A plot of the chromatic zeros for a long finite strip of this type is shown in Fig. [2^. From our exact 
analytic results we find that the locus B includes arcs and a short interval on the real axis 2.078191 < q < q c 
where q c = 2.099012. 
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Figure 23: Chromatic zeros for the width L y = 4 strip of the honeycomb lattice with free boundary conditions, (FBC y , FBC X ) 
and length m = 14 (i.e., n = 118). 

11 Numerical Values of W 

It is of interest to use our new exact results to study further the approach of W to the L y — > oo limit, 
i.e., the infinite 2D honeycomb lattice. This extends our previous study in |p^| . In Table || we list values 
of W(hc,L y x oo, BC y , BC X , q) (which, for this range of q, are independent of BC X ), for free and periodic 
(cylindrical) transverse boundary conditions. These are compared with Monte Carlo measurements of W for 
the full 2D square lattice, W(hc, q) from |50| We also list the ratios 

Rw(A,L v ,BC v ,q)= W ^ A B q ^ q) (11.1) 



for the present case, A = he. In []3j| it was proved for a general lattice A that the ratio ( [11. l[ ) approaches 
unity monotonically for FBC y . For PBC y , this approach is nonmonotonic, but more rapid. One sees 
from Table || that for L y = 4 and moderate values of q, say 5 or 6, the agreement of W(hc, L y , q) for the 
infinite-length, finite- width open strips with the respective values W(hc, q) for the infinite honeycomb lattice 
is very good; the ratios are unity to within a few per cent. For the same width L y = 4 and similar q 
values, the agreement for the cylindrical strip is excellent; the ratios are unity to within a few parts in 10 5 . 
The agreement, as measured by these ratios, gets better as q increases. This shows that calculations on 
infinite-length, finite-width strips of the honeycomb lattice with cylindrical boundary conditions provide an 
excellent approximation to the exact function W(hc, q) for the full 2D honeycomb lattice. This is a very 
useful result since the latter function has never been calculated exactly. We also give some further values of 
W or \W\ in Table |. 



37 



Table 3: Values of W(hc,L y X oo, FBC y , BC X , q), denoted W(hc,L y ,F,q) for short, and W(hc,L y X oo, PBC y , BC x ,q), 
denoted W(hc, L y , P,q) for short, with any BCx- These are compared with the values for the full 2D honeycomb lattice, 
W(hc,q) = W {he, oo X oo, q) for 3 < q < 10. For each value of q, the quantities in the upper line are identified at the top and 
the quantities in the lower line are the values of Rw{hc,L y ,BC y ,q). 



q 


W(hc,2,F,q) 


W(hc,3,F,q) 


W(hc,A,F,q) 


W{hc,A,P,q) 


W(hc,q) 


3 


1.82116 
1.097 


1.76567 
1.064 


1.73864 
1.047 


1.67079 
1.0065 


1.6600(5) 
1 


4 


2.79468 
1.073 


2.72942 
1.048 


2.69737 
1.036 


2.60443 
1.0002 


2.6038(7) 
1 


5 


3.78389 
1.057 


3.71448 
1.038 


3.68025 
1.028 


3.57963 
1.0000 


3.5796(10) 
1 


6 


4.77760 
1.046 


4.70568 
1.031 


4.67013 
1.023 


4.56513 
0.9999 


4.5654(15) 
1 


7 


5.77336 
1.039 


5.69973 
1.026 


5.66327 
1.019 


5.55529 
0.9999 


5.5556(17) 
1 


8 


6.77028 
1.034 


6.69539 
1.0225 


6.65826 
1.017 


6.54810 
1.0000 


6.5479(20) 
1 


9 


7.76793 
1.030 


7.69208 
1.020 


7.65443 
1.015 


7.54259 
1.0000 


7.5424(22) 
1 


10 


8.76607 
1.027 


8.68945 
1.018 


8.65140 
1.013 


8.53822 
1.0000 


8.5386(25) 
1 



Table 4: Values of W(hc, L y X oo, BC y , BC X , q) for low integral q and for respective q c , if such a point exists, where BC y and 
BC X denote the transverse and longitudinal boundary conditions. The notation n in the q c column means that there is no q c 
for this family, i.e., B q does not cross the positive real q axis. 



Ly 


BCy 


BC X 


\w q=0 \ 


\W q =i\ 


\W q =2\ 


Qc 


w q=qc 


2 


F 


F 


1.495 


1 


1 


n 




3 


F 


F 


1.70 


1.26 


1 


2 


1 


4 


F 


F 


1.81 


1.39 


1.09 


2.10 


1.06 


2 


F 


(T)P 


1.495 


1.19 


1 


2 


1 


3 


F 


(T)P 


1.70 


1.35 


1 


2 


1 


4 


P 


F 


2.05 


1.56 


1.19 


2.25 


1.12 
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12 Conclusions 



In this paper we have presented exact calculations of the partition function of the g-state Potts model 
and the special case of the T = Potts antiferromagnet (chromatic polynomial) on a variety of strips of the 
honeycomb lattice with fixed width and arbitrarily great length. We have computed the exact free energy and 
discussed the thermodynamics for the ferromagnetic and antiferromagnetic cases, including the singularities 
at the zero-temperature critical point of the fcrromagnet. These results were compared with those obtained 
from exact free energy calculations for infinite-length strips of the square and triangular lattices and the 
differences explained. The W values calculated from the finite-width, infinite-length strips were shown to 
approach the values for the full 2D honeycomb lattice reasonably rapidly. This was especially dramatic 
for the L y = 4 strip with periodic transverse (cylindrical) boundary conditions, and the exact results for 
this strip thus provide a very accurate approximation to the W function for the full 2D honeycomb lattice. 
This is useful since there is no exact solution for W(q) on the honeycomb lattice for q > 2. Generalizing 
q and temperature to complex values and taking the limit of infinite length, we have studied the locus of 
singularities B of the free energy in the q and u planes, and for the T = antiferromagnet, the locus B q . 



Acknowledgment: The research of R. S. was supported in part by the U. S. NSF grant PHY-97-22101. 

13 Appendix 

13.1 A Remark on Generating Functions 

In (3^, |3^, [37| the methodology on generating function for chromatic polynomials for recursive families of strip 
graphs was given. For generality, the open strip graphs were taken to be of the symbolic form G m = IH m , 
i.e., m-fold repetitions of a subgraph H connected to an initial end graph /. It was proved in |35j that the 
denominator of the generating function is independent of the end graph /. In |36| , a study was carried out 
of inhomogeneous strip graphs where I ^ H. However, if / = H, one can simplify the general formalism 
slightly by redefining the label m via a shift by one unit. Thus, if I = H, then the strip denoted symbolically 



as IH m is identical to H m+1 . The generating function is given in general by equations of the form (9.1) 



fl9.4|) . Now in the convention (which is denoted here as CI) of |35|, |36|, |37j], the zcroth-order term, Aq^o 
in the expansion ( |9.1[ ) is the chromatic polynomial for the subgraph /. If / = H , as in the homogeneous 
strips considered here, then it is simpler to use a convention, which we denote C2, according to which the 
zeroth-order term in the expansion is the chromatic polynomial for one end (say the left-hand end) of the 
strip. For example, for the width L y strip of the square lattice, in conventions CI and C2, A sq x y ,o would 
be the chromatic polynomial for a column of squares of height L y and the line graph Tl , respectively, i.e., 
A sq ,L y ,o = q(q — 1)(<7 2 — 3g + 3) iH_1 and A aqi L q = <z(<7 — l) Ly_1 - To distinguish these two conventions, let 
us, in this subsection, append a prime to the numerator defined with the latter convention C2, i.e., 

Af(G s ,q,z)' = Y, A V.A<1)*' ■ ( 13 - L1 ) 

Thus, for m > 1, the m'th term in the expansion of the generating function F(G S , q, z) 1 is the same as the 
(m — l)'th term in the expansion of the generating function T(G S , q, z). While the denominators are the 
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same, the coefficient functions in the numerators are related according to 

A' Gs!j = Ag.j-1 + X; ■,>'<; ., for j > 1 . (13.1.2) 

Similar relations hold for the generating functions of the full Potts model partition function and cquivalcntly 
the Tutte polynomial. 

We give some simple illustrations next. For the L y — 2 strip of the square lattice with free boundary 
conditions and CI, 

q(g - 1)£> 4 

TP ^ 2 ~ 1-D 4 z (13 - L3) 



while with convention C2 

g(g-i 

1 - D iZ 

For the L y = 3 strip of the square lattice with free boundary conditions and labelling convention CI feq 



= ■ (13.1.4) 



Ap, sg3 , = q(q - l)Dj (13.1.5) 

A P , sq3>1 = -q(q - 1)V - Gq 2 + 13-? - 11) (13.1.6) 

A'p sq3 , Q = q(q - 1) 2 (13.1.7) 

4>.«a,i = -?(9 - 1)(? 2 - 3« + 1) (13.1.8) 

bp sq3 ,i = ~(q - 2)(q 2 -3q + 5) (13.1.9) 

b P , sq z,2 = {q~ 1)(<Z 3 ~ 6q 2 + 13« - 11) . (13.1.10) 

Similarly, with convention CI, the generating functions for the Tutte polynomials for the L y = 2 open 
strip of the square and triangular lattices are 



while for convention C2, 



and in both cases 



(y + x + x 2 +x 3 )- yx 3 z ,„.„> 

7>9. 2 — -, -, — j— ; j 2\ i 2~2 (^id.l.llj 

1 — (y + t + x + x z )z + yx z 



and 

x(x + l) 2 + 2xy + y(y + 1) — x 3 y 2 z 
T ' tn ' 2 = l-[{x + l) 2 + y(y + 2)}z + x 2 y 2 z 2 
while with convention C2, they have the simpler forms 



(13.1.12) 



r' - + y(1 ~ x)z ri3 i 131 

L T,sq.2 — -. i 2\ i 2 2 ^lO.l.lOJ 

H 1 — (y + t + x + x )z + yx z z* 



and 

x - y{yx -y-l)z 
1 - [(ar + l) 2 + y{y + 2)}z + x 2 y 2 z 2 



L T,tri,2 — -, \f„ | ^ 2 i „./„. i oM~ i „.2„.2_2 (XH.X.Xt) 



and so forth for other lattice strips. We shall use convention C2 here. 
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13.2 Connection Between Potts Model Partition Function and Tutte Polyno- 
mial 

The formulas relating the Potts model partition function Z(G, q, v) and the Tutte polynomial T(G, x, y) were 
given in |l4j and hence we shall be brief here. The Tutte polynomial of G, T(G, x, y), is given by 

T(G, x,y)=J2( x ~ l) fe(G,) - fe(G) (y - 1) C(G ' } (13.2.1) 

G'CG 

where k(G'), e(G'), and n(G') = n(G) denote the number of components, edges, and vertices of G' , and 

c(G') = e(G') + fc(G') - n(G') (13.2.2) 
is the number of independent circuits in G 1 . For the graphs of interest here, k(G) = 1. Now let 

x = l + - (13.2.3) 

v 

and 

y = a = v + l (13.2.4) 

so that 

g=(x-l)(y-l). (13.2.5) 

Then 

Z(G, q, v) = (x- l) k ^ (y - l) n ^T(G, x, y) . (13.2.6) 
For a planar graph G the Tutte polynomial satisfies the duality relation 

T(G,x,y)=T(G*,y,x) (13.2.7) 

where G* is the (planar) dual to G. As discussed in |l4j , the Tutte polynomial for recursively defined graphs 
comprised of m repetitions of some subgraph has the form 

T(G m ,x,y) =5> T)G>j (A T , GJ ) m . (13.2.8) 

3=1 

13.3 Open L y = 2 Strip of the Honeycomb Lattice 

The generating function for the Tutte polynomial for the open strip of the honeycomb lattice comprised of 
m bricks, denoted S m , is 

oo 

r T (S m ,x,y,z) = Y / T (S m ,x,y)z m . (13.3.1) 



m=0 



We calculate 



This yields 



where 



r , c , x + y(l-x)z ., 
T T {S, x, y, z) = -— — 5- — — , , n — ; ~ • (13.3.2) 

1 — (x^ + x A + x z + x + y + l)z + yx^z z 



X 



T,S,(1,2) 



2 



1 J • y + l± VR^s 



x + x + x + x - 



(13.3.3) 



R TS = x 8 + 2x 7 + 3x 6 - 2yx A + 4x 5 + 2yx 3 + 5x 4 + 2yx 2 

+4a; 3 + y 2 + 2yx + 3x 2 + 2y + 2x + 1 . (13.3.4) 
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13.4 Cyclic and Mobius L y = 2 Strips of the Honeycomb Lattice 

We calculate 

6 

T(L ro ,x,y)=^ CT , iJ (A T)L , i ) m 



and 



where 



T(ML m , X,y) = ^2 C T, ML,j(^T,ML,jY 
3=1 



^T,ML,j — \r,L,j , .7 = 1) — )6 



At,l,i = 1 

At,L,2 = £ 2 



A 



T,L,(3,4) 



2.x + y + 2 ± V^T34 



with 



i? T34 = x 4 + 4e 3 - 2x 2 y + 8a; 2 + 4xy + 8x + y 2 + Ay + 4 

At,L,5 = At,S,1 , At,L,6 = At,S,2 

where At,s,i an d Xt,s,2 were given above. 

It is convenient to extract a common factor from the coefficients: 



ct.g, 



C T,G, 3 

x-1 



G = L, ML 



(13.4.1) 

(13.4.2) 

(13.4.3) 
(13.4.4) 
(13.4.5) 

(13.4.6) 

(13.4.7) 
(13.4.8) 

(13.4.9) 



Of course, although the individual terms contributing to the Tutte polynomial are thus rational functions of 
x rather than polynomials in x, the full Tutte polynomial is a polynomial in both x and y. We have 

c T ,L,j = c z ,L,j for 1 < j < 6 . (13.4.10) 

Thus, in terms of the variables x and y, e.g., &r L,j = q— l = xy~x — y for j = 2, 3, 4, etc. Further, 

ct,mlj = cz,ML,j for 1 < j < 6 . (13.4.11) 

13.5 Special Values of Tutte Polynomials for Strips of the Honeycomb Lattice 

For a given graph G = (V,E), at certain special values of the arguments x and y, the Tutte polynomial 
T(G,x,y) yields quantities of basic graph-theoretic interest 0)-|O|. We recall some definitions: a spanning 
subgraph was defined at the beginning of the paper; a tree is a connected graph with no cycles; a forest is 
a graph containing one or more trees; and a spanning tree is a spanning subgraph that is a tree. We recall 
that the graphs G that we consider are connected. Then the number of spanning trees of G, Nst(G), is 



Nst(G) = T(G, 1, 1) , 
the number of spanning forests of G, Nsf(G), is 

Nsf(G) = T(G, 2, 1) , 



(13.5.1) 



(13.5.2) 
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the number of connected spanning subgraphs of G, Ncssg(G), is 

N C ssg(G) =T(G, 1,2) , 
and the number of spanning subgraphs of G, Nssg(G), is 

Nssg(G) =T(G,2,2) . 

An elementary theorem (e.g., jl7|]) is that 

N SSG (G) = 2 e < G > . 



(13.5.3) 



(13.5.4) 



(13.5.5) 



From our calculations of Tutte polynomials, we find that for the L y = 2 strip of the honeycomb lattice 
with free boundary conditions, S m , 



NsT(S m ) = 



1 



4V2 



(3 + 2V2 ) m+1 - (3 - 2V2 r +1 



N SF (S m ) 



(31 + 8V15 )[4(4 + V15 )} m - (31 - 8V15 )[4(4 - Vl5 



Ncssg(S 



7 + V4T\ m+1 /7-V41 



m+l 



NsSG(Sra) = 2 e(5 '" } = 2 5m+1 



(13.5.6) 

(13.5.7) 

(13.5.8) 
(13.5.9) 



For the L y = 2 cyclic and Mobius strips of the honeycomb lattice, L m and ML m , we find, with r\Q = ±1 
for L m and ML m , respectively: 



N ST (Gm) = m 



-2 VG + (3 + 2V2 ) m + (3 - 2V2 ) r 



for G m = L m , ML m 



(13.5.10) 



N S F(G m ) = r, G (l~2 2m )- 



11 + V105 Y /ll - V105 



-{[4(4 + Vl5 )] m + [4(4 - Vl5 )] m }for G m = L m , ML m (13.5.11) 



NcssG(L m ) = N CSS G(ML m ) -Am -I 



= -2(m+l) + 



1 + m 1 



1 + m{ 1 - 
1 



41/ x - 2 
7 - \/4l\ m 



^VssG(im) = NssG(ML m ) = 2< L ^ = 2 5m . 



(13.5.12) 
(13.5.13) 
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Since T(G m ,x,y) grows exponentially as m — ► oo for the families G m = S m and L m for (x,y) = (1, 1), 
(2,1), (1,2), and (2,2), it is natural to define the corresponding constants 

z set ({G})= lim 71(G)- 1 In N set (G) , set = ST, SF, CSSG, SSG (13.5.14) 

n(G) — >oo 

where, as above, the symbol {G} denotes the limit of the graph family G as n(G) — > oo (and the z here 
should not be confused with the auxiliary expansion variable in the generating function (13.3.1) or the Potts 
partition function Z(G,q,v).)) General inequalities for these were given in ]I3| |. Our results yield 

zst({G}) = i hi(3 + 2V2 ) ~ 0.440687 for G = S,L,ML (13.5.15) 
z SF {{G}) = - In U(4 + VT5 )] - 0.862433 for G = S,L,ML (13.5.16) 



4 

zcssg{{G}) = \ - 0-475585 for G = S,L,ML (13.5.17) 

and 

zssg{{G}) = (5/4) In 2 ~ 0.866434 for G = S,L,ML (13.5.18) 

Another comparison of interest is the ratio of z$t for these L y = 2 strips with z$t for the full 2D honeycomb 
lattice, which has the value J6(J 

1 3 /3 

z ST (hc) = -z ST (tri) = (1 - 5~ 2 + 7~ 2 - ll" 2 + 13~ 2 - ...) = 0.807664868... (13.5.19) 

2 2tt 

namely, 

ZST ^} ~ 0.545631 . (13.5.20) 

ZST{hc) 

In the case of zst, it is also of interest to compare the results with a general upper bound [ |io[ 

^st({G}) < zsT,u PP er({G}) = In A eff ({G}) (13.5.21) 



where the effective coordination number (degree) was defined above in eq. (4.4). For this comparison we 
define the ratio 

r ST ({G}) = ZST[{G ^ rx , ■ (13.5.22) 

Z ST, upper yX^S) 

We find 

r ST({G}) = ln( , 3 1 + . 2 5 ^ • fOT G = S,L, ML (13.5.23) 

4 Ml) 

We compare these results with others that we have obtained for infinite-length strips of other lattices |l4], [l(| 
in Table ||. We recall that the strip of the triangular lattice is constructed by starting with a strip of the 
square lattice and adding edges joining the lower left to upper right vertex of each square. For uniformity, 
when comparing the values of zst({G}) to upper bounds, we use the general bound ( 13.5.21 ). It should be 
remarked, however, that in contrast to the strip graphs of the honeycomb lattice considered here, the cyclic 
L y = 2 strip graphs of the square and triangular lattice and the infinite-length open strips of these lattices 
are A-regular, and hence one can apply a somewhat more restrictive upper bound, as discussed in j67j]. 

The (planar) dual of L m can be described as follows. In mathematical graph theory, the "join" G + H 
of the graphs G and H is defined as the graph with the vertex set V = Vq + Vh and edge set E comprised 
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Table 5: z s ({G}) for the L x — > oo limit of the width L y = 2 strips of the (i) honeycomb lattice, (ii) square lattice, and (iii) 
triangular lattice. 



z s {{G}) 


{G} = he 


{G} = sq 


{G} = tri 


zst{{G}) 
r ST {{G}) 


(1/4) ln(3 + 2^2 ) 
= 0.440687 
0.4809 


(l/2)ln(2 + V3 ) 
= 0.658479 
0.5994 


(l/2)ln[(7 + 3V5 )/2] 
= 0.962424 
0.6942 


zsf{{G}) 


(1/4) ln[4(4 + ^15 )] 
= 0.862433 


(l/2)ln[2(2 + V3 )] 
= 1.005053 


(l/2)ln[2(3 + 2V2 )] 
= 1.227947 


zcssg{{G}) 


(1/4)M(7 + V41 )/2] 
= 0.475585 


(l/2)ln[(5 + Vl7)/2] 
= 0.758832 


(l/2)ln[2(3 + 2V2 )] 
= 1.227947 


zssg({G}) 


(5/4) In 2 
= 0.866434 


(3/2) In 2 
= 1.039721 


21n2 
= 1.386294 



of the edges of G and H together with edges joining each vertex of G to each vertex of H. Further, the 
complete graph K p is defined as the graph with p vertices such that each vertex is connected by an edge to 
every other vertex. The complement K p of K p is the graph composed of p disjoint vertices and no edges. 
Now consider the join K 2 + C m , where C m is the circuit graph with m vertices. Change each of the edges 
connecting the two vertices of the K 2 to the vertices of the circuit graph C m to double edges. This yields a 
multigraph that we shall denote as K2(+) e 2C m , where the symbol G(+) e ^H means that one replaces each 
edge joining a vertex of G with a vertex of H by i such edges. Then the (planar) dual of L* n is 

L* m = K 2 (+) e2 C m . (13.5.24) 

The planar dual of S m is the graph composed of the join of a single vertex K\ with the line graph T m 
modified so that each edge connecting the external vertex to the (m — 2) internal vertices of the line graph 
is replaced by a 4-fold multiple edge, and the edges connecting this external vertex to the two end vertices 
of the line graph are replaced by 5-fold multiple edges. Thus, our calculations of various graph-theoretic 



quantities for L m and S m also determine those for L* m and by eq. (13.2.7). Note that this relation does 
not apply for the Mobius strips since they are not planar. 

References 

[1] R. B. Potts, Proc. Camb. Phil. Soc. 48 (1952) 106. 
[2] F. Y. Wu, Rev. Mod. Phys. 54 (1982) 235. 
[3] G. D. Birkhoff, Ann. of Math. 14 (1912) 42. 
[4] H. Whitney, Ann. of Math. 33 (1932) 688. 

[5] P. W. Kasteleyn and C. M. Fortuin, J. Phys. Soc. Jpn. 26 (1969) (Suppl.) 11. 

[6] C. M. Fortuin and P. W. Kasteleyn, Physica 57 (1972) 536. 

[7] W. T. Tutte, Can. J. Math. 6 (1954) 80. 

[8] W. T. Tutte, J. Combin. Theory 2 (1967) 301. 



45 



[9] W. T. Tutte, "Chromials", in Lecture Notes in Math. v. 411 (1974) 243; Graph Theory, vol. 21 of 
Encyclopedia of Mathematics and Applications (Addison- Wesley, Menlo Park, 1984). 



[10 
[11 

[12 
[13 

[14 

[is; 

[16 

[17; 

[18 
[19 
[20 
[21 

[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 
[33 
[34 



N. L. Biggs, Algebraic Graph Theory (2nd ed., Cambridge Univ. Press, Cambridge, 1993). 

D. J. A. Welsh, Complexity: Knots, Colourings, and Counting, London Math. Soc. Lect. Note Ser. 186 
(Cambridge University Press, Cambridge, 1993). 

B. Bollobas, Modern Graph Theory (Springer, New York, 1998). 

R. Shrock, in the Proceedings of the 1999 British Combinatorial Conference, BCC99, Discrete Math., 



in press (cond-mat/9908387) 



R. Shrock, Physica A 283 (2000) 388. 

H. Kluepfel and R. Shrock, YITP-99-32; H. Kluepfel, Stony Brook thesis (July, 1999). 



S.-C. Chang and R. Shrock, Physica A, in press ( cond-mat/0004181 ) 



S.-C. Chang and R. Shrock, |cond-mat/0007505 



M. Aizenman and E. H. Lieb, J. Stat. Phys. 24 (1981) 279. 
Y. Chow and F. Y. Wu, Phys. Rev. B36 (1987) 285. 
R. C. Read, J. Combin. Theory 4 (1968) 52. 

R. C. Read and W. T. Tutte, "Chromatic Polynomials", in Selected Topics in Graph Theory, 3, eds. L. 
W. Beineke and R. J. Wilson (Academic Press, New York, 1988.). 

R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 5165. 

N. L. Biggs, R. M. Damerell, and D. A. Sands, J. Combin. Theory B 12 (1972) 123. 
N. L. Biggs and G. H. Meredith, J. Combin. Theory B20 (1976) 5. 
N. L. Biggs, Bull. London Math. Soc. 9 (1976) 54. 

S. Beraha, J. Kahane, and N. Weiss, J. Combin. Theory B 27 (1979) 1. 

S. Beraha, J. Kahane, and N. Weiss, J. Combin. Theory B 28 (1980) 52. 

R. C. Read, in Proc. 3rd Caribbean Conf. on Combin. and Computing (1981). 

R. C. Read, Proc. 5th Caribbean Conf. on Combin. and Computing (1988). 

R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 6791. 

R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 1342, 2733, 4111. 

R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 3935. 

J. Salas and A. Sokal, J. Stat. Phys. 86 (1997) 551. 



R. Shrock and S.-H. Tsai, Phys. Rev. E58 (1998) 4332; |cond-mat/9808057 

46 



[35] M. Rocek, R. Shrock, and S.-H. Tsai, Physica A252 (1998) 505. 

[36] M. Rocek, R. Shrock, and S.-H. Tsai, Physica A259 (1998) 367. 

[37] R. Shrock and S.-H. Tsai, Physica A259 (1998) 315. 

[38] R. Shrock and S.-H. Tsai, J. Phys. A Lett. 32 (1999) L195. 

[39] R. Shrock and S.-H. Tsai, Phys. Rev. E60 (1999) 3512. 

[40] R. Shrock and S.-H. Tsai, Physica A 275 (1999) 429. 

[41] R. Shrock and S.-H. Tsai, J. Phys. A 32 (1999) 5053. 

[42] N. L. Biggs, LSE report LSE-CDAM-99-03 (May 1999), to appear. 

[43] R. Shrock, Phys. Lett. A261 (1999) 57. 

[44] N. L. Biggs and R. Shrock, J. Phys. A (Letts.) 32 (1999) L489. 

[45] S.-C. Chang and R. Shrock, Stony Brook preprints YITP-SB-99-50,58. 



[46] S.-C. Chang and R. Shrock, Physica A, in press (cond-mat/0004161) 



[47] S.-C. Chang and R. Shrock, Physica A, in press ( cond-mat /0005232 ) 



[48] S.-C. Chang and R. Shrock, |cond-mat/0007491 . 

[49] V. Matveev and R. Shrock, J. Phys. A 28 (1995) 1557. 

[50] V. Matveev and R. Shrock, J. Phys. A 29 (1996) 803. 

[51] V. Matveev and R. Shrock, Phys. Lett. A204 (1995) 353. 

[52] D. Kim and R. Joseph, J. Phys. C 7 (1974) L167. 

[53] T. W. Burkhardt and B. W. Southern, J. Phys. A 11 (1978) L247. 

[54] R. Shrock and S.-H. Tsai, J. Phys. A 30 (1997) 495. 

[55] J. Salas, J. Phys. A 31 (1998) 5969. 

[56] H. Feldmann, R. Shrock, and S.-H. Tsai, J. Phys. A 30 (1997) L663. 

[57] H. Feldmann, R. Shrock, and S.-H. Tsai, Phys. Rev. E57 (1998) 1335. 

[58] M. E. Fisher, Lectures in Theoretical Physics (Univ. of Colorado Press, 1965), vol. 7C, p. 1. 

[59] R. Abe, Prog. Theor. Phys. 38 (1967) 322. 

[60] C. N. Yang and T. D. Lee, Phys. Rev. 87 (1952) 404. 

[61] T. D. Lee and C. N. Yang, ibid. 87 (1952) 410. 

[62] H. Feldmann, A. J. Guttmann, I. Jensen, R. Shrock, and S.-H. Tsai), J. Phys. A 31 (1998) 2287. 



47 



[63] C. N. Chen, C. K. Hu, and F. Y. Wu, Phys. Rev. Lett. 76 (1996) 169. 

[64] F. Y. Wu, G. Rollet, H. Y. Huang, J.-M. Maillard, C. K. Hu, and C. N. Chen, Phys. Rev. Lett. 76 
(1996) 173. 

[65] V. Matveev and R. Shrock, Phys. Rev. Phys. Rev. E54 (1996) 6174. 

[66] F. Y. Wu, J. Phys. A 10 (1977) L113. 

[67] R. Shrock and F. Y. Wu, J. Phys. A 33 (2000) 3881. 



48 



